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I review the origin and properties of electromagnetic fields produced in heavy ion collisions. 
The field strength immediately after a collision is proportional to the collision energy and 
reaches ~ at RHIC and <~ lOmJ at LHC. I demonstrate by explicit analytical calculation 
that after dropping by about one-two orders of magnitude during the first fm/c of plasma 
expansion, it freezes out and lasts for as long as quark-gluon plasma lives as a consequence 
of finite electrical conductivity of the plasma. Magnetic field breaks spherical symmetry in 
the direction perpendicular to the reaction plane and therefore all kinetic coefficients are 
anisotropic. I examine viscosity of QGP and show that magnetic field induces azimuthal 
anisotropy on plasma flow even in spherically symmetric geometry. Very strong electromag- 
netic field has an important impact on particle production. I discuss the problem of energy 
loss and polarization of fast fermions due to synchrotron radiation, consider photon decay 
induced by magnetic field, elucidate J/ip dissociation via Lorentz ionization mechanism and 
examine electromagnetic radiation by plasma. I conclude that all processes in QGP are 
affected by strong electromagnetic field and call for experimental investigation. 
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1. ORIGIN AND PROPERTIES OF ELECTROMAGNETIC FIELD 
§1. Origin of magnetic field 

We can understand the origin of magnetic field in heavy-ion collisions by considering collision 
of two ions of radius R with electric charge Ze (e is the magnitude of electron charge) at impact 
parameter b. According to the Biot and Savart law they create magnetic field that in the center- 
of-mass frame has magnitude 

B~ 1 Ze^ (1.1) 

and points in the direction perpendicular to the reaction plane (span by the momenta of ions) . Here 
7 = \/snn /2m7v is the Lorentz factor. At RHIC heavy-ions are collided at 200 GeV per nucleon, 
hence 7 = 100. Using Z = 79 for Gold and b ~ R A « 7 fm we estimate eB « ml ~ 10 18 G. To 
appreciate how strong is this field, compare it with the following numbers: the strongest magnetic 
field created on Earth in a form of electromagnetic shock wave is ~ 10 7 G [2], magnetic field of a 
neutron star is estimated to be 10 10 — 10 13 G, that of a magnetar up to 10 15 G [1]. It is perhaps 
the strongest magnetic field that have ever existed in nature. 

It has been known for a long time that classical electrodynamics breaks down at the critical 
(Schwinger) field strength F = m^/e [3]. In cgs units the corresponding mag netic field is 10 13 G. 
Because m^/rrie = 280, electromagnetic fields created at RHIC and LHC are well above the critical 
value. This offers a unique opportunity to study the super-strong electromagnetic fields in labora- 
tory. The main challenge is to identify experimental observables that are sensitive to such fields. 
The problem is that nearly all observables studied in heavy-ion collisions are strongly affected both 
by the strong color forces acting in quark-gluon plasma (QGP) and by electromagnetic fields often 
producing qualitatively similar effects. An outstanding experimental problem thus is to separate 
the two effects. In sections 2-7 I examine several processes strongly affected by intense magnetic 
fields and discuss their phenomenological significance. But first, in this section, let me derive a 
quantitative estimate of electromagnetic field. 

Throughout this article, the heavy-ion collision axis is denoted by z. Average magnetic field 
then points in the y-direction, see Fig. 1 and Fig. 7. Plane xz is the reaction plane and b is the 
impact parameter. 



FIG. 1: Heavy- ion collision geometry as seen along the collision axis z. Adapted from [7]. 



§2. Magnetic field in vacuum 



A. Time- dependence 



To obtain a quantitative estimate of magnetic field we need to take into account a realistic 
distribution of protons in a nucleus. This has been first done in [4]*. Magnetic field at point r 
created by two heavy ions moving in the positive or negative z-direction can be calculated using 
the Lienard-Wiechert potentials as follows 



with R a = r — r a (t), where sums run over all Z protons in each nucleus, their positions and 
velocities being r a and v a . The magnitude of velocity v a is determined by the collision energy 



approximation, assuming that protons travel on straight lines before and after the scattering. This 
is a good approximation since baryon stopping is a small effect at high energies. Positions of 
protons in heavy-ions can be determined by one of the standard models of the nuclear charge 
density p(r a ). Ref.[4] employed the "hard sphere" model, while [7] used a bit more realistic Wood- 
Saxon distribution. 

Numerical integration in (1.3) including small contribution from baryon stopping yields for 
magnetic field the result shown in Fig. 2 as a function of the proper time r = (t 2 — z 2 ) x l 2 . Evidently, 
magnetic field rapidly decreases as a power of time, so that after first 3 fm it drops by more than 




(1.2) 



(1.3) 



yJsNN and the proton mass m p , v 2 = 



1 — (2m p / ^/snn) . These formulas are derived in the eikonal 



* In the case of high-energy pp collisions, magnetic field was first estimated in [5] who also pointed out a possibility 
of formation of W-condensate [5, 6]. 
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three orders of magnitude 

io E 




(a) 



FIG. 2: Magnetic field B = By (multiplied by e) at the origin r — produced in collision of two gold ions 
at beam energies (a) y/s^N = 62 GeV and (b) v /ijv r iv = 200 GeV. Adapted from [4]. Note, that eB is the 
same in Gauss and Lorcntz-Hcavisidc units in contrast to B. 



B. Event-by- event fluctuations in proton positions 

Nuclear charge density p provides only event-averaged distribution of protons. The actual 
distribution in a given event is different form p implying that in a single event there is not only 
magnetic field along the y-direction, but also other components of electric and magnetic fields. 
This leads to event-by-event fluctuations of electromagnetic field [7]. Shown in Fig. 3 are electric 
and magnetic field components at t = at the origin (denoted by a black dot in Fig. 1) in AuAu 
collisions at ^/snn = 200 GeV. 

Fig. 3 clearly shows that although on average the only non- vanish component of the field is B y , 
which is also clear from the symmetry considerations, other components are finite in each event 
and are of the same order of magnitude 

(\B X \) » (\E X \) » (\E y \) . (1.4) 

To appreciate the magnitude of electric field produced in heavy-ion collisions note that E ~ = 
10 21 V/cm. The corresponding intensity is 10 39 W/cm 2 which is instructive to compare with the 
power generated by the most powerful state-of-the-art lasers: 10 23 W/cm 2 . 
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FIG. 3: The mean absolute value of (a) magnetic field and (b) electric field at t = and r = as a function 
of impact parameter b for AuAu collision at ^snn — 200 GeV. 



Electromagnetic fields created in heavy-ion collisions were also examined in more elaborated 
approaches in [8-10]. They yielded qualitatively similar results on electromagnetic field strength 
and its relaxation time. 



§3. Magnetic field in quark-gluon plasma 

A. Lienard-Wiechert potentials in static medium 

In the previous section I discussed electromagnetic field in vacuum. A more realistic estimate 
must include medium effects. Indeed, state-of-the art phenomenology of quark-gluon plasma (QGP) 
indicates that strongly interacting medium is formed at as early as 0.5 fm/c. Even before this time, 
strongly interacting medium exist in a form of Glasma [11, 12]. Therefore, a calculation of magnetic 
field must involve response of medium determined by its electrical conductivity. It has been found 
in the lattice calculations that the gluon contribution to electrical conductivity of static quark-gluon 
plasma is [14] 

T 

a = (5.8 ± 2.9) — MeV , (1.5) 

J c 

where T is plasma temperature and T c its critical temperature. This agrees with [15], but is 
at odds with an earlier calculation [16]. It is not clear whether (1.5) adequately describes the 
electromagnetic response of realistic quark-gluon plasma because it neglects quark contribution 
and assumes that medium is static. Theoretical calculations are of little help at the temperatures 
of interest since the perturbation theory is not applicable. In absence of a sensible alternative I 



8 



will use (1.5) as a best estimate of electrical conductivity. If medium is static then T is constant 
as a function of time t. The static case is considered in this section, while in the next section I 
consider expanding medium. 

In medium, magnetic field created by a charge e moving in z-direction with velocity v is a 
solution of the following equations 

V-B = 0, VxE = ~, (1.6) 

dE 

V • E = eS(z - vt)S(b) , VxB = — + aE + evzS(z — vt)8(b) , (1.7) 

where we used the Ohm's law j = aE to describe currents induced in the medium. Position of the 
observation point is specified by the longitudinal and transverse coordinates z and b, r = zz + b. 
Taking curl of the second equation in (1.7) and substituting eqs. (1.6) we get 

V 2 B=°^ r +a°^- evV x [zd(z - vt)5(b)} . (1.8) 

The particular solution reads 

f't /"CO /* 

B(z, b,t)= / dt' / dz d 2 bG(z- z',b-b',t- t') evV' x [z5(z' - vt')5(b')] , (1.9) 

JO J-oo J 

where the Green's function G(z — z',b — b', t — t') satisfies the following equation 

V 2 G - ^ - a— = S(z - z')S(b - b')S(t - t>) , (1.10) 

which is solved by 

G(z- z',b-b',t-t') 

d2k± e i(b-b')-k ± f°° ^z e ik z (z-z') j°° du_ e -iu(t-t') 1 ^ ^ L1 ^ 



(2vr) 2 J_ 00 2ir J_ OQ 2ir k z + k\ — W — iuja 

where k = k z z + k±. Plugging this into (1.9) and substituting for the expression in the square 
brackets in (1.9) its Fourier image we obtain 

d 2 k± j b . fc f 00 dk z Ak _ /Z f 00 duj __ iujt ik x z 



B(z, b, t) = 2irev / ^^e lb ' fc± / -^ e ^ 2 / — e~ tuJt = 5(u - k z v) , 

^ ' J (2vr) 2 7_ 00 2vr 2ir k 2 + k\ - u 2 - iua V ; 

(1.12) 

d 2 k± jb.k ± f°° iwt i0JZ / v ik± x z 

e / o_ e e . .9 /„.9 , 7.2 .9 ■ ._■ \ L - LA ) 



(2tt) 2 J_ oq 2tt uj 2 jv 2 + k\ — ui 2 — iua 

We are interested in the y-component of the field. Noting that (k± x z) ■ y = —k± cos</>, where 4> 
is the azimuthal angle in the transverse plane, and integrating over d 2 k± we derive 

eB v = - / s^idfsH^e"^')^, (1.14) 

J-oo 
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where we introduced notation 



— 1 ~'\/^- e M- ( L15 ) 



e(u) is the dielectric constant of the plasma with the following frequency dependence 



io~ 



e(w) = l + — . (1.16) 



UJ 



Eq. (1.14) is actually valid for any functional form of e(uj) [13], which can be easily verified by 
using electric displacement D instead of E in eqs. (1.7). In this case (1.16) can be viewed as a low 
frequency expansion of e(uj). Magnetic field in this approximation is quasi-static. Therefore, we 
could have neglected the second time derivative in (1.8) and then keeping only the leading powers 
of uj we would have derived (1.14) with s 2 = iujo. After integration over uj this gives (1.21). Let us 
take notice of the fact that neglecting the second time derivative in (1.8) yields diffusion equation 
for magnetic field in plasma. 

It is instructive to compare time-dependence of magnetic field created by moving charges in 
vacuum and in plasma. In vacuum, setting a = in (1.13) and integrating first over uj and then 
over k± gives 

eB = ^a cm(62 + ^ 7 _ z)2)3/2 , (1.17) 

where we used v ~ 1. This coincides with (1.3) for a single proton when we take R a = b + (z — vt)z. 
Consider field strength (1.17) at the origin z = 0. At times t < 6/7 the field is constant, while at 
t S> 6/7 it decreases as oc 1/t 3 . At the time t « b the ratio between these two is 

-Bo 1 , 
/ = -,«!, (1.18) 

which is a very small number (~ 1(T 6 at RHIC). 

In matter a > 0. Let me write the modified Bessel function appearing in (1.14) as follows 

r , , [°° J 1 (xb)x 2 dx 

sKi(sb) = j Q n x2 ^ g2 . (1.19) 

Substituting (1.19) into (1.14) and using (1.16) we have (v = 1) 



eB = ^y dx dw ? ' e™*-** (1.20) 

V JO J -oc X 2 - IUJG 



' - J x {xb)x 2 M ,_ t) 

Closing the contour in the lower half-plane of complex uj picks a pole at uj = —ix 2 /a. We have 

&B = 2a^^ /■« dxx 2 J 1 {xb)e- x ^ = y p^e'W^ . (1.21) 
a J 2(t — z) 
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At z = this function vanishes at t = and t — > oo and has maximum at the time instant t = b 2 a/8 
which is much larger than 6/7. The value of the magnetic field at this time is 

32e~ 2 a PTn 



b 3 a 



(1.22) 



(Here e is the base of natural logarithm). This is smaller that the maximum field in vacuum 



Br, 



32e 



-2 



;i.23) 



Bq airy 

but is still a huge field. We compare the two solutions (1.17) and (1.21) in Fig. 4. We see that in 
a conducting medium magnetic field stays for a long time. 

One essential component is still missing in our arguments - time-dependence of plasma prop- 
erties due to its expansion. Let us now turn to this problem. 

ZeB,fm~ 2 



0.001 




t.fm 



FIG. 4: Relaxation of magnetic field at z = in vacuum (blue), in static conducting medium at a = 5.8 MeV 
(red) and at a — 16 MeV (brown) and in the expanding medium (green). Units of B is fm~ 2 w 2mJ. 
6 = 7 fin, Z = 79 (Gold nucleus), 7 = 100 (RHIC). 



B. Magnetic field in expanding medium 

So far I treated quark-gluon plasma as a static medium. In expanding medium temperature 
and hence conductivity are functions of time. In Bjorken scenario [17], expansion is isentropic, i.e. 
nV = const, where n is the particle number density and V is plasma volume. Since n ~ T 3 and 
at early times expansion is one- dimensional V ~ t it follows that T oc i -1 / 3 . (Eventually, we will 
consider the midrapidity region z = 0, therefore distinction between the proper time and t is not 
essential). Eq. (1.5) implies that a ~ t^ 1 / 3 . I will parameterize conductivity as follows 

"( ( >="»(^) 1/3 ' (L24) 
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where I took to ~ 0.5 fm to be the initial time (or longitudinal size) of plasma evolution. Suppose 
that plasma lives for 10 fm/c and then undergoes phase transition to hadronic gas at T c . Then 
employing (1.5) we estimate cfq ~ 16 MeV. Let me define another parameter that I will need in 
the forthcoming calculation: 

" = £»«™ 

Magnetic field in expanding medium is still governed by (1.8). As was explained in the preceding 
subsection, time-evolution of magnetic field is quasi-static, which allows me to neglect the second 
time derivative. Let me introduce a new "time" variable p as follows 

p = (1 + t/t f 3 - 1 • (1.26) 

Field B(z,b,p) satisfies equation 

V 2 B = 0^- - evV x {z5[z - vt(p)]S(b)} , (1.27) 

where 

t(p) = t [(p + l) 3/4 -l]. (1.28) 

Its solution can be written as 

B(z,b,p)= f f 'dpi ( dz I d 2 b<5(z- z',b-b',p- p')evV x {zS[z - vt(p)]5(b)} , (1.29) 

JO J-oo J 

in terms of the Green's function (5(z — z',b — b', p — p') satisfying 

V 2 © - = -5(z - z')8(b - b')5(p - p') , (1.30) 

To solve this equation we represent <3 as three-dimensional Fourier integral with respect to the 
space coordinates and Laplace transform with respect to the "time" coordinate: 

.(,. M = / <^<— > I ^e"^ , (1.31, 

with the contour C running parallel to the imaginary axis to the right of all integrand singularities. 
Now I would like to write the expression in the curly brackets in (1.29) also as Fourier-Laplace 
expansion. To this end we calculate 

d 2 b / dz dpe- i( - k ^ b+k ^e- xp 5(z-vt(p))S(b) (1.32) 

J-oo Jo 

poo 

= / dpe -ik.«to[(P+i) s/4 -i] e -Ap. (1.33) 
Jo 
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Therefore, 

evV x {z5[z - vt{p)]5(b)} = ev j e i(k±-b+k z z) J g A P - fe± x £ ^ _ ( L34 j 
Substituting (1.31) and (1-34) into (1.29) we obtain upon integration over the volume and time 

B(2 , M ./ / c A , A w , (1 ,5, 

where 6* is the step-function. Taking consequent integrals over A and k z gives 

B(0,M = ^ / ^%^' 5 r — *± >< * rVe-'^^ 1 ^- 11 e"^ 1 ^^ (1.36) 

= Z \ 7^%e^ b ^ik ± x z / d p'e-^ p - p) 4=e ^ • (1-37) 

(3 J (2vr) 2 2vr 7 v 7 ^ 7 

Consider now Integrating over azimuthal angle (j) and then over k± as in (1.13), (1.14) yields 
eB y (0,b,p) = ^^ d(C 5/2 e~e * , (1.38) 

where Q = p — p' . 

The results of a numerical calculation of (1.38) are shown in Fig. 4. We see that expansion of 
plasma tends to increase the relaxation time, although this effect is rather modest. We conclude 
that due to finite electrical conductivity of QGP, magnetic field essentially freezes in the plasma 
for as long as plasma exists. Similar phenomenon, known as skin-effect, exists in good conduc- 
tors placed in time-varying magnetic field: conductors expel time dependent magnetic fields form 
conductor volume confining them into a thin layer of width 5 ~ u^ 1 / 2 on the surface. 



C. Diffusion of magnetic field in QGP 

The dynamics of magnetic field relaxation in conducting plasma plasma can be understood in 
a simple model [18]. Suppose at some initial time t = magnetic field B(0, r) permits the plasma. 
The problem is to find the time-dependence of the field at t > 0. In this model, the field sources 
turn off at t = and do not at all contribute to the field at t > 0. Electromagnetic field is governed 
by the following equations 

j = aE (1.39) 
V-S = 0, (1.40) 



VxB = j, 
„ „ dB 
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that lead to the diffusion equation for B, after we neglect the second time derivative as discussed 
before 

V* B = ff §. (1.41) 

For simplicity we treat electrical conductivity a as constant. Initial condition at t = reads 

B(0,r) = B e~^ , (1.42) 

where the Gaussian profile is chosen for illustration purposes and R is the nuclear radius. Solution 
to the problem (1.41), (1.42) is 

B(t,r) = J dV'B(0,r')G(t,r-r'). (1.43) 

where the Green's function is 

1 T r-2 1 

(1.44) 



G(t, r) = . exp 



{Ant/a) 3 / 2 1 I 4t/<7_ ' 
Integrating over the entire volume we derive 



2 



B(t, r) = B exp 



b 



2 



(1.45) 



R 2 + 4t/a ' I R 2 +At/a 

It follows from (1.45) that as long as t <C trelax, where t re i ax is a characteristic time t rc i ax = R 2 & /4 
magnetic field B is approximately time-independent. This estimate is the same as the one we 
arrived at after (1.21). 

In summary, magnetic field in quark-gluon plasma appears to be extremely strong and slowly 
varying function of time for most of the plasma life-time. At RHIC it decreases from eB « (2.5m 7r ) 2 
right after the collision to eB ps at t « 5 fm, see Fig. 4. This has a profound impact on 

all the processes occurring in QGP. 

D. Schwinger mechansim 

Schwinger mechanism of pair production [3] is operative if electric field exceeds the critical 
value of m 2 /e, where m is mass of lightest electrically charged particle. Indeed, in order to excite 
a fermion out of the Dirac sea, electric force eE must do work along the path di satisfying 

/ eE-d£'>2m. (1.46) 
Jo 
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If E = const, then E > m/ie. The maximal value of I is the fermion Compton's wavelength 
I ~ A = 1/m implying that the minimum (or critical) value of electric field is 

Ec = ^. (1.47) 

Notice that in stronger fields I ~ m/eE < A. Fig. 3 indicates that electron-positron pairs are 
certainly produced at RHIC. An important question then is the role of these pairs in the elec- 
tromagnetic field relaxation in plasma. There are two associated effects: (i) before e~e + pairs 
thermalize, they contribute to the Foucault currents, (ii) after they thermalize, their density con- 
tributes to the polarization of plasma in electric field, and hence to its conductivity. 

Since space dimensions of QGP are much less than A e = 380 fm, it may seem inevitable that 
space-dependence of electric field (in addition to its time-dependence) has a significant impact 
on the Schwinger process in heavy-ion collisions. However, his conclusion is premature. Indeed, 
suppose that electric field is a slow function of coordinates. Then E(r) ~ E(0) + r ■ Vi2(0). Work 
done by electric field is 

[ eE(0)-d£' + f (r • V)e.E(0) • dl' ~ eEl + e ^°) ; (i. 48 ) 
Jo Jo A 

where A is length scale describing space variation of electric field. In order that contribution of space 

variation to work be negligible, the second term in the r.h.s. of (1.48) must satisfy e£ 2 E(0)/X <C m. 

Employing the estimate i ~ m/eE(0) that we obtained after (1.47) implies m/eE(0)X <C 1. 

Following [19] I define the inhomogeneity parameter 

7=^ (1.49) 

that describes the effect of spatial variation of electric field on the pair production rate. For 
electrons m = 0.5 MeV in QGP A ~ 0.5 fm at eE ~ m% we have 7 ~ 0.01. Therefore, somewhat 
counter-intuitively , electric field can be considered as spatially homogeneous. The same conclusion 
can be derived from results of [20]. Schwinger mechanism in spatially-dependent electric fields was 
also discussed in [21, 22]. 

In view of smallness of 7 one can employ the extensive literature on Schwinger effect in time- 
dependent spatially-homogeneous electric fields. As far as heavy-ion physics is concerned, the most 
comprehensive study has been done in [23-25] who developed an approach to include the effect of 
back reaction. They argued that time-evolution of electric field can be studied in adiabatic approx- 
imation and used the kinetic approach to study the time-evolution. Their results are exhibited in 
Fig. 5. Similar results were obtained in [26]. We observe that response time of the current density 
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of Schwinger pairs ~ 10 4 fm/c is much larger than the plasma lifetime ~ 10 fm/c and therefore no 
sizable electric current is generated. 




FIG. 5: Time-dependence of electric field due to the Schwinger mechanism back-reaction and the corre- 
sponding electric current density of Schwinger pairs. Dimensionless time variable is defined as r = t/X. For 
electrons A = 380 fm. Plasma undergoes phase transition at about t = 1/38. Adapted from [23]. 

In summary, strong electric field is generated in heavy ion collisions in every event, but 
averages to zero in a large event ensemble. This field exceeds the critical value for electrons and 
light u, d quarks. However, during the plasma lifetime no significant current of Schwinger pairs is 
generated. 
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2. FLOW OF QUARK-GLUON PLASMA IN STRONG MAGNETIC FIELD 

§1. Azimuthal asymmetry 

Magnetic field is known to have a profound influence on kinetic properties of plasmas. Once the 
spherical symmetry is broken, distribution of particles in plasma is only axially symmetric with 
respect to the magnetic field direction. This symmetry however is not manifest in the plane span 
by magnetic field and the impact parameter vectors, viz. xy-plane in Fig. 1. Charged particles 
moving along the magnetic field direction y are not influenced by the magnetic Lorentz force 
while those moving the xz-plane (i.e. the reaction plane) are affected the most. The result is 
azimuthally anisotropic flow of expanding plasma in the xy-plane even when initial plasma geometry 
is completely spherically symmetric. The effect of weak magnetic field on quark-gluon plasma 
flow was first considered in [27] who argued that magnetic field is able to enhance the azimuthal 
anisotropy of produced particles up to 30%. This conclusion was reached by utilizing a solution of 
the magneto-hydrodynamic equations in weak magnetic field. 

A characteristic feature of the viscous pressure tensor in magnetic field is its azimuthal 
anisotropy. This anisotropy is the result of suppression of the momentum transfer in QGP in 
the direction perpendicular to the magnetic field. Its macroscopic manifestation is decrease of 
the viscous pressure tensor components in the plane perpendicular to the magnetic field, which 
coincides with the reaction plane in the heavy-ion phenomenology. Since Lorentz force vanishes 
in the direction parallel to the field, viscosity along that direction is not affected at all. In fact, 
the viscous pressure tensor component in the reaction plane is twice as small as the one in the 
field direction. As the result, transverse flow of QGP develops azimuthal anisotropy in presence of 
the magnetic field. Clearly, this anisotropy is completely different from the one generated by the 
anisotropic pressure gradients and exists even if the later are absent. In fact, because spherical 
symmetry in magnetic field is broken, viscous effects in plasma cannot be described by only two 
parameters: shear r\ and bulk viscosity (. Rather the viscous pressure tensor of magnetoactive 
plasma is characterized by seven viscosity coefficients, among which five are shear viscosities and 
two are bulk ones. 
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§2. Viscous pressure in strong magnetic field 

A. Viscosities from kinetic equation 

Generally, calculation of the viscosities requires knowledge of the strong interaction dynamics 
of the QGP components. However, in strong magnetic field these interactions can be considered 
as a perturbation and viscosities can be analytically calculated using the kinetic equation [28, 29]. 
To apply this approach to QGP in strong magnetic field we start with kinetic equation for the 
distribution function / of a quark flavor of charge ze is 

p^f = zeB^^u u + C[f,...] (2.1) 

where C is the collision integral and B^ v is the electro-magnetic tensor, which contains only mag- 
netic field components in the laboratory frame. Ellipsis in the argument of C indicates the dis- 
tribution functions of other quark flavors and gluons (I will omit them below). The equilibrium 
distribution reads 

where U(x) is the macroscopic velocity of the fluid, p^ = mu^ is particle momentum, /3 = 1/T and 
p is the mass density. Since oc u^, the first term on the r.h.s. of (2.1) as well as the collision 
integral vanishe in equilibrium. Therefore, we can write the kinetic equation as an equation for Sf 

V„ = zeB^^Uv + C[6f] (2.3) 

where 5f is a deviation from equilibrium. Differentiating (2.2) we find 

d„fo = -fo^P%U x (x) (2.4) 

Since U x = (7y,7yV) and p x = (e,p) = {~f v m,^ v mv) it follows that 



Thus, in the comoving frame 



Substituting (2.6) in (2.3) yields 



dufo\v=o = fo^Pud u V (2.6) 



-^P v V U v = zeB^^u„ + C[5f] (2.7) 
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where I defined 

V liV = ±(d l M, + d v V li ) (2.8) 

and used u^u v d^V u = u^u u V^ v . 

Since the time-derivative of /o is irrelevant for the calculation of the viscosity I will drop it from 
the kinetic equation. All indices thus become the usual three- vector ones. To avoid confusion we will 
label them by the Greek letters from the beginning of the alphabet. Introducing b a p = B~ 1 e a p~ f B 1 
we cast (2.7) in the form 

^P a u fi V aP f = -zeBb a ^^-- £ - C[Sf] . (2.9) 
The viscous pressure generated by a deviation from equilibrium is given by the tensor 

-IW = J PaPpSf^ (2.10) 

Effectively it can be parameterized in terms of the viscosity coefficients as follows (we neglect the 
bulk viscosities) 



na/3 = E^ F i? ( 2 - n ) 
n=0 



(n) 

where the linearly independent tensors are given by 



y(0) 

Oi/3 


= (3b a bp - S a p) (b^sV^s - • V^j 


(2.12a) 


v m 

oi/3 


= 2V aj3 + SapVysb-ybs - 2V aj b 1 bp - 2Vp 1 b 1 b a + (b a bp - 5 a p)V ■ 


V + b a bpV l5 b 1 b$.l2b) 


V ap 


= 2{V ai bp 1 + Vpjb ai - V^sba^bpbg) 


(2.12c) 


V a/3 


= Vaybp^ + Vp^ba-y — V y s b a s 6 Q7 bp bs — V^sbp^babs 


(2.12d) 




= 2(Vjsbasb ai bpbs + Vysbp^b a bs) . 


(2.12e) 



For the calculation of the shear viscosities rj n , n = 1, . . . , 4 we can set V • V = and V a pb a bp = 0. 

(n) 

Let us expand of to the second order in velocities in terms of the tensors as follows 

4 
n=0 

Then, substituting (2.13) into (2.11) and requiring consistency of (2.10) and (2.11) yields 

Vn = ~hJ £vignd * p (2 - 14) 

This gives the viscosities in the magnetic field in terms of the deviation of the distribution function 
from equilibrium. Transition to the non-relativistic limit in (2.14) is achieved by the replacement 
e ->• m [28]. 
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B. Collisionless plasma 

In strong magnetic field we can determine g n by the method of consecutive approximations. 
Writing 6f = 8fW + 6f<® and substituting into (2.9) we find 

^V a ,f = -zeBb a ^ mi l^ f{2)) l + C[5f^ . (2.15) 

Here I assumed that the deviation from equilibrium due to the strong magnetic field is much larger 
than due to the particle collisions. The explicit form of C is determined by the strong interaction 
dynamics, but drops off the equation in the leading oder. The first correction to the equilibrium 
distribution obeys the equation 

1 <9<5/W 1 

T^PaVpVapfo = -zeBb a/3 Vf3 — - . (2.16) 



Using (2.13) we get 

dv a 



bapvp 8 ^ = 2b aP vp 9n Vffv^ . (2.17) 



n=0 

Substituting (2.17) into (2.16) and using (2.12) yields: 

e 



TzeB 



PaVpVapfo = -2bp v v a v v [g l {2V a p-2Vf3 1 b~ i b a ) + 2g 2 Vp~ i b~ 1 b a 

+93{V ay b^ + V Pl b ai - Vy S b a b 5 ) + 2g A V lS b^b a bs)} , (2.18) 



where I used the following identities b a pb a = b a pbp = b a pv a vi3 = 0. Clearly, (2.18) is satisfied only 
if 9i = 92 = 0. Concerning the other two coefficients, we use the identities 

bapb^ = b 1 b a -5 ai b 2 , (2.19a) 
e a/ 3 7 £5ec = S a s(Sp £ S^ — S^S ie ) — d ae (6^6^ — 6/3^6^)+ Sa^i^SpsSje — S^Sjs) (2.19b) 

that we substitute into (2.18) to derive 

~Wz~eB pav ^ fo = 9&V a pb a bp - W a pv a bp(b ■ v)} + 2g 4 V a pv a b f} {b ■ v) . (2.20) 
Since p a = ev a we obtain 

9s = ^ = —^-. (2.21) 

y 2 ATzeB v ' 

Using (2.2), (2.21) in (2.14) in the comoving frame (of course ry n 's do not depend on the frame 
choice) and integrating using 3.547.9 of [51] we derive [30] 

% = ^Mj^. (2 . 22) 

13 K 2 (f3m)2zeB [ } 
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The non-relativistic limit corresponds to m > T in which case we get 




pT 



(2.23) 



2zeB 



In the opposite ultra-relativistic case m <C T (high-temperature plasma) 




2nT 2 



(2.24) 



zeB 



where n = p/m is the number density. 



C. Contribution of collisions 



In the relaxation-time approximation we can write the collision integral as 



C[6f] = -u6f, 



(2.25) 



where v is an effective collision rate. Strong field limit means that 



LOB S> V , 



(2.26) 



where cob = zeB/e is the synchrotron frequency. Whether v itself is function of the field depends 
on the relation between the Larmor radius r# = vt/cob, where vt is the particle velocity in the 
plane orthogonal to B and the Debye radius ro- If 



then the effect of the field on the collision rate v can be neglected [28]. Assuming that (2.27) is 
satisfied the collision rate reads 



where at is the transport cross section, which is a function of the saturation momentum Q s [31, 32]. 
We estimate at ~ a 2 /Q 2 , with Q s ~ 1 GeV and n = P/T with pressure a 2 s P ~ 1 GeV/fm 3 we get 
v ~ 40 MeV. Inequality (2.26) is well satisfied since eB ~ m 2 [4, 8] and m is in the range between 
the current and the constituent quark masses. On the other hand, applicability of the condition 
(2.27) is marginal and is very sensitive to the interaction details. In this section we assume that 
(2.27) holds in order to obtain the analytic solution. Additionally, the general condition for the 
applicability of the hydrodynamic approach £ = \ jv ^ L, where t is the mean free path and L is 
the plasma size is assumed to hold. Altogether we have t\d <C tb <C i <C L. 



tb > r D 



(2.27) 



v = nvat , 



(2.28) 
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Equation for the second correction to the equilibrium distribution 5f^ follows from (2.15) after 
substitution (2.25) 



Now, plugging 



<V (1) = faV^+giVfflvaVp, (2.30a) 

*/3 + 92V a/ 3 



¥ 2) = [giV$+g 2 vM]v a vp, (2.30b) 



into (2.29) yields 
2zeB 



{giftVpaba^vpv^ - 2Vp a b ai vpv 1 {v ■ b)\ + 2g 2 Vp a b Q . 1 vpv 1 {v ■ b)} 

= -vgz {-ZVpaba-fVpV-f - GVpabayVpV^V • b)} , (2.31) 



where I used 54 = 2g%. It follows that 



91 = J = ^- (2 - 32) 



With the help of (2.28), (2.2), (2. 14) we obtain [30] 



= V1 8 pig t T^ K 7/2 (Pm) 
Vl 4 5^ {zeBfm 1 / 2 K 2 {pm) ' { ' ' 



§3. Azimuthal asymmetry of transverse flow: a simple model 

To illustrate the effect of the magnetic field on the viscous flow of the electrically charged 
component of the quark-gluon plasma I will assume that the flow is non-relativistic and use the 
Navie-Stokes equations that read 



dt dxp J dx a dxp ' 
where H a p is the viscous pressure tensor, p = mn is mass-density and P is pressure. I will 
additionally assume that the flow is non-turbulent and that the plasma is non-compressible. The 
former assumption amounts to dropping the terms non-linear in velocity, while the later implies 
vanishing divergence of velocity 

V-V = 0, (2.35) 

Because of the approximate boost invariance of the heavy-ion collisions, we can restrict our atten- 
tion to the two dimensional flow in the xy plane corresponding to the central rapidity region. 



22 



The viscous pressure tensor in vanishing magnetic field is isotropic in the xy-plane and is given 

by 

^,-v(^ + ^)-J r - V -\. (2.36, 



\dx/3 dx c 



Vy X Vyy 



where the superscript indicates absence of the magnetic field. In the opposite case of very strong 
magnetic field the viscous pressure tensor has a different form (2.11). Neglecting all r\ n with n > 1 
we can write 



I -Vyy 




| = 2% | 


2 y xx 










V 


2Vyy 




V 


Vyy 



V% = vo\ = 2?7o , (2.37) 

\ 2V yy ) \ Vyy) 

where we also used (2.35). Notice that H"^. = jll^ = \li xx indicating that the plasma flows in 
the direction perpendicular to the magnetic field with twice as small viscosity as in the direction 
of the field. The later is not affected by the field at all, because the Lorentz force vanishes in the 
field direction. Substituting (2.37) into (2.34) we derive the following two equations characterizing 
the plasma velocity in the strong magnetic field [30] 

dV x dP d 2 V x dV v dP d 2 V v . , 

Additionally, we need to set the initial conditions 

V x \ t=0 = <pi(x,y), V y \ t=0 = <p 2 (x,y) . (2.39) 
The solution to the the problem (2. 38), (2. 39) is 

V x (x,y,t)= dx'ip 1 {x',y)Gi(x-x',t) - - f dt' H dx'Gi {x -x',t- t') dP ^^ Q ; 

J-oo 2 P JO J-oo 2 V x 

(2.40a) 

V y (x,y,t)= [°° dy'<p 2 (x, 1 /)G 1 (y-y',t)-± f dt' f°° dy'G 1 (y - y' ,t - t >) dP M ^) . 

J-oo P Jo J-oc Off 

(2.40b) 



Here the Green's function is given by 



G k {y,t) = -2==e-^i (2.41) 
V4:ira 2 kt 



and the diffusion coefficient by 
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Suppose that the pressure is isotropic, i.e. it depends on the coordinates x,y only via the radial 
coordinate r = \J x 2 + y 2 ; accordingly we pass from the integration variables x' and y' to r' in 
(2.40a) and (2.40b) correspondingly. At later times we can expand the Green's function (2.41) 
in inverse powers of t. The first terms in the r.h.s. of (2.40a) and (2.40b) are subleading and we 
obtain 

1 f* , f°° , 1 dP(r,t-s) 

V x (x,y,t)Bs — / ds dr 

PJo J-oo V2ira 2 s dr 

= -- f ds^L= [P(R(s), t-s)- P(0, t-s)}, (2.43a) 
P Jo V27ra / s 

and by the same token 

V y (x,y,t)^-- [ ds—L=[P(R(s),t-s)-P(0,t-s)], (2.43b) 
P Jo y/Aita^s 

where R(t) denotes the boundary beyond which the density of the plasma is below the critical 
value. We observe that V x /V y = y/2. Consequently, the azimuthal anisotropy of the hydrodynamic 
flow is [30] 

V 2 - V 2 1 - i 1 

V 2 + V 2 -l + \ "3- (2 ' 44j 

Since I assumed that the initial conditions and the pressure are isotropic, the azimuthal asymmetry 
(2.44) is generated exclusively by the magnetic field. 

— 1/2 

We see that at later times after the heavy-ion collision, flow velocity is proportional to rj , 
where 770 is the finite shear viscosity coefficient, see (2.40a) and (2.40b). If the system is such that 
in absence of the magnetic field it were azimuthally symmetric, then the magnetic field induces 
azimuthal asymmetry of 1/3, see (2.44). The effect of the magnetic field on flow is strong and 
must be taken into account in phenomenological applications. Neglect of the contribution by the 
magnetic field leads to underestimation of the phenomenological value of viscosity extracted from 
the data [33-35]. In other words, the more viscous QGP in magnetic field produces the same 
azimuthal anisotropy as a less viscous QGP in vacuum. 

A model that I considered in this section to illustrate the effect of the magnetic field on the 
azimuthal anisotropy of a viscous fluid flow does not take into account many important features 
of a realistic heavy-ion collision. To be sure, a comprehensive approach must involve numerical 
solution of the relativistic magnetohydrodynamic equations with a realistic geometry. A potentially 
important effect that I have not considered here, is plasma instabilities [37, 38], which warrant 
further investigation. 
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The structure of the viscous stress tensor in very strong magnetic field (2.37) is general, model 
independent. However, as explained, the precise amount of the azimuthal anisotropy that it gen- 
erates cannot be determined without taking into account many important effects. Even so, I draw 
the reader's attention to the fact that analysis of [27] using quite different arguments arrived at 
similar conclusion. Although a more quantitive numerical calculation is certainly required before a 
final conclusion can be made, it looks very plausible that the QGP viscosity is significantly higher 
than the presently accepted value extracted without taking into account the magnetic field effect 
[33-35] and is perhaps closer to the value calculated using the perturbative theory [39, 40]. 
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3. ENERGY LOSS AND POLARIZATION DUE TO SYNCHROTRON RADIATION 
§1. Radiation of fast quark in magnetic field 

General problem of charged fermion radiation in external magnetic field was solved in [41-43]. 
It has important applications in collider physics, see e.g. [44, 45]. In heavy-ion phenomenology, 
synchrotron radiation provides one of the mechanisms of energy loss in quark-gluon plasma, which 
is an important probe of QGP [46, 47] J 

A typical diagram contributing to the synchrotron radiation, i.e. radiation in external magnetic 
field, by a quark is shown in Fig. 6 [18]. This diagram is proportional to (eB) n , where n is the 
number of external field lines. In strong field, powers of eB must be summed up, which may be 
accomplished by exactly solving the Dirac equation for the relativistic fermion and then calculating 
the matrix element for the transition q — > q + g. Such calculation has been done in QED for some 
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FIG. 6: A typical diagram contributing to the synchrotron radiation by a quark. 

special cases including the homogeneous constant field and can be readily generalized for gluon 
radiation. Intensity of the radiation can be expressed via the invariant parameter x defined as 



X 2 = -^p- { F, u ff = ^p-( P x B) 2 (3.1) 



where the initial quark 4-momentum is p 11 = (e,p), Z q is the quark charge in units of the absolute 
value of the electron charge e. At high energies, 

Z„Be 

X«lf— ■ 3.2 
B c m 

The regime of weak fields corresponds to x ^ 1> while in strong fields j( > 1. In our case, 
eB/eB c « {ra^jra,^ 2 3> 1 (at RHIC) and therefore x>f In terms of x, spectrum of radiated 



' Synchrotron radiation in chromo-magnetic fields was discussed in [48-50]. 
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gluons of frequency oj can be written as [41] 



= -a s C F r ^|^ O °Ai(O^+0 + ^^ 1/2 )Ai'(x)} , (3.3) 



where / is the intensity, 

and e' is the quark's energy in the final state. Ai is the Ayri function. Eq. (3.3) is valid under the 
assumption that the initial quark remains ultra-relativistic, which implies that the energy loss due 
to the synchrotron radiation Ae should be small compared to the quark energy itself Ae <C e. 
Energy loss by a relativistic quark per unit length is given by [44] 

de f°° dl m 2 X 2 r 4 + 5xx 3 / 2 + 4 X 2 x 3 

dl=-J d ^ = asCF — J, (l + x *3/ 2 )4 Ai(x)xdx. (3.4) 
In two interesting limits, energy loss behaves quite differently. At rj = ip = we have [44] 

de _ 2a s hC F (Z q eB) 2 e 2 

dl ~ 3rr? ' X<<1 ' (3 ' 5a) 

J = -0.37 ash-^Cp {Z q eBe) 2/s , (3.5b) 

In the strong field limit energy loss is independent of the quark mass, whereas in the weak field 
case it decreases as m~ 4 . Since % oc H, limit of x *C 1 corresponds to the classical energy loss. 

To apply this result to heavy-ion collisions we need to write down the invariant x in a suitable 
kinematic variables. The geometry of a heavy-ion collision is depicted in Fig. 7. Magnetic field 
B is orthogonal to the reaction plane span by the impact parameter vector b and the collision 
axis (z-axis). For a quark of momentum p we define the polar angle 6 with respect to the z- 
axis and azimuthal angle <p with respect to the reaction plane. In this notation, B = B y and 
p = p z z+p±(x cos if + y sin (p) , where p± = \p\ sin 6 ps esin6*. Thus, (Bxp) 2 = B 2 (p 2 z +p 2 L cos 2 tp). 
Conventionally, one expresses the longitudinal momentum and energy using the rapidity rj as 
e = m± cosh?? and p z = m± sinhr/, where m\ = m 2 + p\. We have 



X 2 = p\ (sinh 2 n + cos 2 p) . (3.6) 



In Fig. 8 a numerical calculation of the energy loss per unit length in a constant magnetic field 
using (3.4) and (3.6) is shown [18]. We see that energy loss of a u quark with p± = 10 GeV is about 
0.2 GeV/fm at RHIC and 1.3 GeV/fm at LHC. This corresponds to the loss of 10% and 65% of 
its initial energy after traveling 5 fm at RHIC and LHC respectfully. Therefore, energy loss due to 
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FIG. 7: Geometry of a heavy-ion collision, p denotes momentum of a fast quark. Two orthogonal planes 
are the reaction plane span by the initial heavy- ion momenta, and the mid-rapidity plane 9 — ir/2, which is 
labeled as 77 = 0. 

the synchrotron radiation at LHC gives a phenomenologically important contribution to the total 
energy loss. 

Energy loss due to the synchrotron radiation has a very non-trivial azimuthal angle and rapidity 
dependence that comes from the corresponding dependence of the x-parameter (3.6). As can be 
seen in Fig. 8(c), energy loss has a minimum at <p = ir/2 that corresponds to quark's transverse 
momentum p± being parallel (or anti-parallel) to the field direction. It has a maximum at 92 = 0, it 
when p± is perpendicular to the field direction and thus lying in the reaction plane. It is obvious 
from (3.6) that at mid-rapidity rj = the azimuthal angle dependence is much stronger pronounced 
than in the forward/backward direction. Let me emphasize, that the energy loss (3.4) divided by 
m 2 , i.e. de/(dlm 2 ) scales with %. In turn, \ is a function of magnetic field, quark mass, rapidity 
and azimuthal angle. Therefore, all the features seen in Fig. 8 follow from this scaling behavior. 



28 




> 

QJ 

a 



CO 



3.5 
3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



eB=15/77~,?/=^=0 



10 



20 
Pr(GeV) 



30 



40 



(a) 



(b) 



u,eB 




. >7=1 


* 


Q 


17=0 \ 

\ Reaction plane 


\ 


v 


/ 

✓ 









1 



<] I CO 
I 



0.025 



0.020 



0.015 




(c) 



(d) 



FIG. 8: Energy loss per um/j length by quarks moving in constant external magnetic field as a function of (a) 
Pt at RHIC at 77 = = 0, (b) px at LHC at r\ = 93 = 0, (c) azimuthal angle ip with respect to the reaction 
plane for p T = 5 GeV and i] = 0, 1. (d) Fractional energy loss vs rapidity i] at = 5 GeV and 1/5 = 0, 7r/4. 

§2. Azimuthal asymmetry of gluon spectrum 

In magnetic field gluon spectrum is azimuthally asymmetric. It is customary to describe this 
asymmetry by Fourier coefficients of intensity defined as 

1 



I (if) =111 + 2 v n cos(ncp) 



n=l 



I(ip) cos(ncp) dip . 



(3.7) 



2tt/ 7_ T1 

Azimuthal averaged intensity is I = J_ I(<p)d(p/2ir. In strong fields x ^ 1 an d we can write 

= 0.37 a s C F m 2 x 2/S = 0.37 a s C F ( eJ Bp ± ) 2/3 (sinh 2 7/ + cos 2 ip) 1/3 . (3.8) 

We have 

f 77 I(u?) cos(nup)dir) f n fsinh 2 r? + cos 2 ip) 1 / 3 cos(nio)dio 

(3.9) 



JZ n I(<p)dip /^(sinh 2 77 + cos 2 c^) 1 / 3 ^ 

At 77 = the Fourier coefficients w n can be calculated analytically using formula 3.631.9 of [51] 



n /4 1 n 4 _ n\ 



, if n 6 even ; -u n = , if n G odd . 



(3.10) 
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where B is the Euler's Beta-function. The corresponding numerically values of the lowest harmonics 
are v 2 = 0.25, v 4 = -0.071, v 6 = 0.036, v 8 = -0.022, v 10 = 0.015. 



§3. Polarization of light quarks 

Synchrotron radiation leads to polarization of electrically charged fermions; this is known as 
the Sokolov-Ternov effect [43]. It was applied to heavy-ion collisions in [18]. Unlike energy loss 
that I discussed so far, this is a purely quantum effect. It arises because the probability of the 
spin-flip transition depends on the orientation of the quark spin with respect to the direction of 
the magnetic field and on the sign of fermion's electric charge. The spin-flip probability per unit 
time reads [43] 

W= ^—^U { e ) [l-- 9 (C-v)-—si g n(e q )«;. z )j, (3.11) 

where £ is a unit axial vector that coincides with the direction of quark spin in its rest frame, 
v = p/e is the initial fermion velocity. 

The nature of this spin-flip transition is transparent in the non-relativistic case, where it is 
induced by the interaction Hamiltonian [45] 

« = = < 312 » 

It is seen, that negatively charged quarks and anti-quarks (e.g. u and d) tend to align against the 
field, while the positively charged ones (e.g. u and d) align parallel to the field. 

Let %(4.) be the number of fermions or anti-fermions with given momentum and spin direction 
parallel (anti-parallel) to the field produced in a given event. At initial time t = to the spin- 
asymmetry defined as 

vanishes. Eq. (3.11) implies that at later times, a beam of non-polarized fermions develops a finite 
asymmetry given by [43] 

15 

where 



(l-e-'^), (3.14) 



-r^-^(-) 2 (^7^l (3-15) 



5V3ma s C F \eJ \Z q e\v x B 



30 



is the characteristic time over which the maximal possible asymmetry is achieved. This time is 
extremely small on the scale of to. For example, it takes only r = 6 • 1CT 8 fm for a d quark of 
momentum p± = 1 GeV at rj = ip = at RHIC to achieve the maximal asymmetry of A m = = 
92%. Therefore, quarks and anti-quarks are polarized almost instantaneously after being released 
from their wave functions. However, subsequent interaction with QGP and fragmentation wash 
out the polarization of quarks. 

A more sensitive probe are leptons weakly interacting with QGP and not undergoing a fragmen- 
tation process. Thus, their polarization can present a direct experimental evidence for the existence 
and strength of magnetic field. In case of muons we can estimate r by replacing a s Cp — > a cm . For 
muons we get r = 0.3 fm/c, which is still much smaller than magnetic field life-time. Observation of 
such a lepton polarization asymmetry is perhaps the most definitive proof of existence of the strong 
magnetic field at early times after a heavy-ion collision regardless of its later time-dependence. 

In summary, energy loss per unit length for a light quark with px = 15 GeV is about 
0.27 GeV/fm at RHIC and 1.7 GeV/fm at LHC, which is comparable to the losses due to in- 
teraction with the plasma. Thus, the synchrotron radiation alone is able to account for quenching 
of jets at LHC with p± as large as 20 GeV. Synchrotron radiation is definetely one of missing 
pieces in the puzzle of the jet energy loss in heavy ion collisions. Quarks and leptons are expected 
to be strongly polarized in plasma in the direction parallel or anti-parallel to the magnetic field 
depending on the sign of their electric charge. 
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4. PHOTON DECAY 

In this section I consider pair-production by photon in external magnetic field [52], which is 
a cross-channel of synchrotron radiation discussed in the previous section. Specifically, we are 
interested to determine photon decay rate w in the process jB — > ffB, where / stands for a 
charged fermion, as a function of photon's transverse momentum kx, rapidity rj and azimuthal 
angle ip. Origin of these photons in heavy-ion collisions will not be of interest to us in this section. 

Characteristic frequency of a fermion of species a of mass m a and charge z a e (e is the absolute 
value of electron charge) moving in external magnetic field B (in a plane perpendicular to the field 
direction) is 

^ B = ^, (4.1) 

e 

where e is the fermion energy. Here - in the spirit of the adiabatic approximation - B is a slow 
function of time. Calculation of the photon decay probability significantly simplifies if motion of 
electron is quasi-classical, i.e. quantization of fermion motion in the magnetic field can be neglected. 
This condition is fulfilled if Hujb *C e. This implies that 

e»VzeB. (4.2) 

For RHIC it is equivalent to e > m^, for LHC s S> Am^. 

Photon decay rate was calculated in [41] and, using the quasi-classical method, in [53]. It reads 

v a cm z 3 a eB /*°° 2(x 3 / 2 + l/x a ) Ai'(s) 
V rn a x a 7(4/^)2/3 x"/4( x 3/2_ 4/Xa) 3/2> 

where summation is over fermion species and the invariant parameter x is defined as 

2 _ a cm^ 3 /p , V \2 _ Qcm^qfc 3 (J R x2 (a a\ 

a ~ ( ^ ^ ~ ^6 ( fc X B > ' ( 4 - 4 ) 

'"-a " l a 

with the initial photon 4-momentum fc M = (faj, fe). With notations of Fig. 7, k = k z z+kj_(x cos 92+ 
ys'mip), where k± = |fc|sin# = hcosmO. Thus, (B x fc) 2 = B 2 (k 2 + Ar^ cos 2 <£>) . Employing 
hio = /cj_ cosh 7] and fc 2 = sinh r\ we write 

>^a = ^ Za& ^ k±\/smh 2 i] + cos 2 ip . (4.5) 

Plotted in Fig. 9 is the photon decay rate (4.3) for RHIC and LHC. The survival probability 
of photons in magnetic field is P = 1 — wAt, where At is the time spent by a photon in plasma. 
Estimating At = 10 fm we determine that photon survives with probability Prhic ~ 97% at 
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(a) (b) 

FIG. 9: Decay rate of photons moving in reaction plane in magnetic field as a function of transverse 
momentum kx'- (a) at RHIC, (b) at LHC. Broken lines from bottom to top give contributions of 7 — ¥ dd 7 
7 — > uu, 7 — > and 7 — > e + e~ channels. Upper solid line is the total rate. 

RHIC, while only -Plhc ~ 80% at LHC. Such strong depletion can certainly be observed in heavy- 
ion collisions at LHC. 

Azimuthal distribution of the decay rate of photons at LHC is azimuthally asymmetric as can 
be seen in Fig. 10 [52]. The strongest suppression is in the B field direction, i.e. in the direction 
orthogonal to the reaction plane. At 77 > 1 the <p dependence of K a is very weak which is reflected 
in nearly symmetric azimuthal shape of the dashed line in Fig. 10. 




Reaction plane 



FIG. 10: Azimuthal distribution of the decay rate of photons at different rapidities at LHC. Only contribution 
of the 7 — ¥ e + e~ channel is shown. 
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To quantify the azimuthal asymmetry it is customary to expand the decay rate in Fourier series 
with respect to the azimuthal angle. Noting that w is an even function of <p we have 

w(<p) = -wq + w n cos(n<f) , w n = — w((p) cos(mp) d<p . (4-6) 

2 n=l nJ-n 

In strong fields x a 3> 1. For example, for 7 — > /U + /U~ at RHIC at ip = rj = and kx = 1 GeV we 
get x M = 19. Therefore, we can expand the rate (4.3) at large x a as [41] 

3 1/6 5r 2 (f) crne ^3 _ a 

7\l^ 1/3 - ( aiT ,y i 2„ , ™ 2,„U/6 ' Xfl ^ 1 ' ^- U 



2V3 77 rV2r(|) ^ maX V3 (sinh^ + cosV) 1 / 6 
At 77 = the Fourier coefficients w n can be calculated analytically using formula 3.631.9 of [51] 

3 2 1 / 3 A 

w 2 k= , 5 5 7T, ^ 2 fc+i = 0, k = 0,1,2,..., (4.8) 

\6 ' 6 / 

where -B is the Euler's Beta-function and A is defined in (4.7). Substituting these expressions into 
(4.6) we find 



W = -Wn 
2 



1 

wo 1 - 

k 

The first few terms in this expansion read 



(4.9) 



1/2 14 \ 

u; = -wo ( 1 - - cos(20) + — cos(4</>) - ... J , (4.10) 

What is measured experimentally is not the decay rate, but rather the photon spectrum. This 
spectrum is modified by the survival probability P which is obviously azimuthally asymmetric. To 
quantify this asymmetry we write using (4.6) 

(°° \ 1 — P 2 

i+y;«2fccos(2^fc) , v 2k = — , (4.11) 
iTi J P wo 

where P = (1 — wAt)^ = 1 — woAt is the survival probability averaged over the azimuthal angle. 
Since woAt <C 1, as can be seen in Fig. 9, we can estimate using (4.7) and (4.8) 

^ 2w 2k 2u, 2fc/w 5 6 2 / 3 r(§) ^ a cm (eB) 2 / 3 z S a /3 
v 2k ~ w At = At — > . • 4.12 

Wo Wo ITT [Kt) I 

In particular, the "elliptic flow" coefficient is [52] 

^^ ^Iff 1 , (4 ,3, 

For example, at kr = 1 GeV and At ~ 10 fm/c one expects v 2 ^ 2% at RHIC and v 2 — 14% at 
LHC only due to the presence of magnetic field. We see that decay of photons in external magnetic 
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field significantly contributes to the photon asymmetry in heavy-ion collisions along with other 
possible effects. 

In summary, I calculated photon pair-production rate in external magnetic field created in 
off-central heavy-ion collisions. Photon decay leads to depletion of the photon yield by a few 
percent at RHIC and by as much as 20% at the LHC. The decay rate depends on the rapidity 
and azimuthal angle. At mid-rapidity the azimuthal asymmetry of the decay rate translates into 
asymmetric photon yield and contributes to the "elliptic flow" . Let me also quote a known result 
that photons polarized parallel to the field are 3/2 times more likely to decay than those polarized 
transversely [41]. Therefore, polarization of the final photon spectrum perpendicular to the field 
is a signature of existence of strong magnetic field. Finally, photon decay necessarily leads to 
enhancement of dilepton yield. 
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5. QUARKONIUM DISSOCIATION IN MAGNETIC FIELD 
§1. Effects of magnetic field on quarkonium 

Strong magnetic field created in heavy-ion collisions generates a number of remarkable effects on 
quarkonium production some of which I will describe in this section. Magnetic field can be treated 
as static if the distance A over which it significantly varies is much larger than the quarkonium 
radius. If At is magnetic field life-time, then A ~ cAt. For a quarkonium with binding energy Sf, and 
radius a s /eb, the quasi-static approximation applies when ef,X/a s S> 1. Estimating conservatively 
A ~ 2 fm we get for J/tp: Eb\/a s ps 23, which is comfortably large to justify the quasi-static 
approximation, where I assumed that is given by its vacuum value. As temperature T increases 
£{, drops. Temperature dependence of is model dependent, however it is certain that eventually it 
vanishes at some finite temperature To. Therefore, only in the close vicinity of To, i.e. at very small 
binding energies, the quasi-static approximation is not applicable. I thus rely on the quasi-static 
approximation to calculate J/tp dissociation [55, 56]. 

Magnetic field has a three-fold effect on quarkonium: 

1. Lorentz ionization. Consider quarkonium traveling with constant velocity in magnetic field 
in the laboratory frame. Boosting to the quarkonium comoving frame, we find mutually 
orthogonal electric and magnetic fields given by Eqs. (5.1), (5. 2). In the presence of an 
electric field quark and antiquark have a finite probability to tunnel through the potential 
barrier thereby causing quarkonium dissociation. In atomic physics such a process is referred 
to as Lorentz ionization. In the non-relativistic approximation, the tunneling probability is 
of order unity when the electric field E in the comoving frame satisfies eE > m 1 / 2 e^ 2 (for 
weakly bound states), where m is quark mass, see (5.24). This effect causes a significant 
increase in quarkonium dissociation rate; numerical calculation for J/tp is shown in Fig. 13. 

2. Zeeman effect. Energy of a quarkonium state depends on spin S, orbital angular momentum 
L, and total angular momentum J. In a magnetic field these states split; the splitting energy 
in a weak field is AM = ^§^gJz, where J z = — J, — J + 1, . . . , J is projection of the total 
angular momentum on the direction of magnetic field, m is quark mass and g is Lande 
factor depending on J, L and S in a well-known way, see e.g. [57]. For example, J/tp with 
S = 1, L = and J = 1 (g ps 2) splits into three states with J z = ±1,0 and with mass 
difference AM = 0.15 GeV, where we used eBo = 15m 2 . Thus, the Zeeman effect leads to 
the emergence of new quarkonium states in plasma. 
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3. Distortion of the quarkonium potential in magnetic field. This effect arises in higher order 
perturbation theory and becomes important at field strengths of order B ~ 37rm 2 /e 3 [58]. 
This is 37r/a times stronger than the critical Schwinger's field. Therefore, this effect can be 
neglected at the present collider energies. 

Some of the notational definitions used in this section: V and P are velocity and momentum 
of quarkonium in the lab frame; M is its mass; p is the momentum of quark or anti-quark in the 
comoving frame; m is its mass; Bo is the magnetic field in the lab frame, E and B are electric and 
magnetic fields in the comoving frame; jl is the quarkonium Lorentz factor; and 7 is a parameter 
defined in (5.19). I use Gauss units throughout the section; note that expressions eB, eE and eB$ 
are the same in Gauss and Lorentz-Heaviside units. 

§2. Lorentz ionization: physical picture 

In this section I focus on Lorentz ionization, which is an important mechanism of J/ip suppres- 
sion in heavy-ion collisions [55, 56]. Before we proceed to analytical calculations it is worthwhile 
to discuss the physics picture in more detail in two reference frames: the quarkonium proper frame 
and the lab frame. In the quarkonium proper frame the potential energy of, say, antiquark (with 
e < 0) is a sum of its potential energy in the binding potential and its energy in the electric field 
—eEx, where x is the electric field direction, see Fig. 11. Since \e\Ex becomes large and negative 
at large and negative x (far away from the bound state) and because the quarkonium potential 
has finite radius, this region opens up for the motion of the antiquark. Thus there is a quantum 
mechanical probability to tunnel through the potential barrier formed on one side by the vanishing 
quarkonium potential and on the other by increasing absolute value of the antiquark energy in 
electric field. Of course the total energy of the antiquark (not counting its mass) is negative after 
tunneling. However, its kinetic energy grows proportionally to eEx as it goes away. By picking up 
a light quark out of vacuum it can hadronize into a L>-meson. 

If we now go to the reference frame where E = and there is only magnetic field B (we can 
always do so since E < B), then the entire process looks quite different. An energetic quarkonium 
travels in external magnetic field and decays into quark-antiquark pair that can late hadronize into 
.D-mesons. This happens in spite of the fact that J/ip mass is smaller than masses of two D-mesons 
due to additional momentum eA supplied by the magnetic field. Similarly a photon can decay into 
electron-positron pair in external magnetic field. 
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§3. Quarkonium ionization rate 

A. Comoving frame 

Consider a quarkonium traveling with velocity V in constant magnetic field -Bo. Let B and 
E be magnetic and electric fields in the comoving frame, and let subscripts || and _L denote field 
components parallel and perpendicular to V correspondingly. Then, 

E\\=0, E ± = lL VxB , (5.1a) 

r B ° V R (F x B ) x V 

B \\ = — y—> B L = 1L — 2 , (5.1b) 

where 7l = (1 — V 2 )~ l l 2 . Clearly, in the comoving frame B ■ E = 0. If quarkonium travels at 
angle <p with respect to the magnetic field in the laboratory frame, then 



B = B, 



b^cos 2 <Ml- 7 2) + 7 2, E = B o7L Vsm<p. (5.2) 

We choose y and x axes of the comoving frame such that B = By and E = Ex. A convenient gauge 
choice is A = —Bxz and ip = —Ex. For a future reference we also define a useful dimensionless 
parameter p [59] 

E 7 L Fsin0 

B ^cos2 0(l- 7 2) + 7 2 

Note, that (i) < p < 1 because B 2 — E 2 = B 2 > and (ii) when quarkonium moves perpendicu- 
larly to the magnetic field Bq, p = V. 



B. WKB method 

I assume that the force binding q and q into quarkonium as a short-range one i.e. {Mei,) l / 2 R <C 1, 
where and M are binding energy and mass of quarkonium, respectively, and R is the radius of 
the nuclear force given by R « (a s /o) 1 / 2 , where a = 1 GeV/fm is the string tension. For example, 
the binding energy of c and c in «//?/> in vacuum is e& = 0.64 GeV <C M/R 2 = Ma/a s « 3 GeV. 
This approximation is even better at finite temperature on account of Ef, decrease. Regarding J/iJj 
as being bound by a short-range force enables us to calculate the dissociation probability w with 
exponential accuracy w ps e~^, independently of the precise form of the quarkonium wave function. 
This is especially important since solutions of the relativistic two-body problem for quarkonium 
are not readily available. 
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It is natural to study quarkonium ionization in the comoving frame [55]. As explained in the 
Introduction, ionization is quantum tunneling through the potential barrier caused by the electric 
field E. In this subsection I employ the quasi-classical WKB approximation to calculate the 
quarkonium decay probability w. For the gauge choice specified in Sec. 5 §3 A quark energy eo 
(eo < m ) in electromagnetic field can be written as 

£o = \J m 2 + (p — eA) 2 + e<p = J m 2 + (p z + eBx) 2 + p 2 + p 2 — eEx . (5-4) 

In terms of £q, quarkonium binding energy is £5 = ui — Eq. To simplify notations, we will set p x = 0, 
because the quark moves constant momentum along the direction of magnetic field. 

The effective potential U(x) = £o(x) — \J m 2 + p 2 corresponding to (5.4) is plotted in Fig. 11. 
We can see that the tunneling probability is finite only if E > 0. It is largest when B = 0. It has 
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FIG. 11: Effective potential U(x) = J m 2 + (p z + eBx) 2 + p%. — eEx — J m 2 + p 2 + p 2 forp x = 0, p z = m/6, 
B = m (except the blue line where B — 0) . The width of the potential barrier decreases with E and increases 
with B. 1 — eo corresponds to the binding energy in units of m. 

been already noted before in [59-61] that the effect of the magnetic field is to stabilize the bound 
state. In spite of the linearly rising potential (at B > E) tunneling probability is finite as the 
result of rearrangement of the QED vacuum in electric field. 

Ionization probability of quarkonium equals its tunneling probability through the potential 
barrier. The later is given by the transmission coefficient 

w = e - 2 fo 1 V z pldy = e -f, (5.5) 

In the non-relativistic approximation one can also calculate the pre-exponential factor, which ap- 
pears due to the deviation of the quark wave function from the quasi-classical approximation. This 
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is discussed later in Sec. 5 §5 B. We now proceed with the calculation of function /. Since B > E 
Eq. (5.4) can be written as [55] 

p 2 x = -e 2 (B 2 -E 2 )(x-x 1 )(x-x 2 ), (5.6) 

where 

e E - p z B T V(e E - Pz B) 2 - - E 2 )(-e 2 + m 2 + p 2 ) 

Xl ' 2 = e{B 2 - E*) • (5 ' 7) 

Define dimensionless variables eo = Eo/m and q = p z /m. Integration in (2.27) gives: 
/ _ y/-4 + l + q 2 (e E - qB) 



m 2 e(B 2 - E 2 ) 

- (B 2 

<B 2 - E 2 f/ 2 1 ^(eo^ - qB) 2 - (B 2 - E 2 )(e 2 + 1 + q 2 ) 



{eoE _ qB ) 2 - (B 2 - E 2 )(-e 2 + l+q 2 ) x J e E - qB + ^(B 2 - E 2 )(-e 2 + 1 + q 2 ) 



(5.8) 

For different q's w = gives the corresponding ionization probabilities. The largest probability 
corresponds to smallest /, which occurs at momentum q m determined by equation [60] 

df(qn 



ImJ 



0. (5.9) 



dq m 

Using (5.8) and parameter p defined in (5.3) we find [55] 

p(ep ~ pq m ) ln / eoP ~ q m + \/j - P 2 V~4 + 1 = V~~4 + 1 + Qm , g 1Q x 

1-P 2 "I V(to- P q m ) 2 -l + P 2 J" \/r^7 

This is an implicit equation for the extremal momentum q m = q m (eo,p). Substituting q m into 
(5.8) one obtains / = f(eo,p), which by means of (5.5) yields the ionization probability. The 
quasi-classical approximation that we employed in this section is valid inasmuch as f(q m ) 3> 1- 

In order to compare with the results obtained in [60] using the imaginary time method, we can 
re- write Eq. (5.10) in terms of an auxiliary parameter To as 

r 0= V /r ^^- g o + 1 ^ ; (5.11a) 
p(e - pq m ) 

tanhr e - pq m , C11U 
=p . (5.11b) 

to eoP - Qm 

Taking advantage of these equations, Eq. (5.8) can be cast into a more compact form 

fm= ^= =[l-e (eo-g ro p)], (5.12) 
eEyJl - p 2 

where we denoted f m = f(q m ). This agrees with results of Ref. [60]. 
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C. Special case: Crossed fields 

An important limiting case is crossed fields E = B. Since also E ± B, see Sec. 5 §3 A, both 
field invariants vanish. Nevertheless, quarkonium ionization probability is finite [60]. This limit 
is obtained by taking p — > 1 in the equations from the previous section. Employing (5.11a) and 
(5.11b) we get the following condition for extremum 

e 2 -l + 2ql-3e q m = 0, (5.13) 

with the solution 

Qm = \ (3e - • (5.14) 

Substituting into (5.12) produces 

tm -\% t*±l±Mt . (5,5) 
3eE eo - q m 

§4. Non-relativistic approximation 

A very useful approximation of the relativistic formulas derived in the previous section is the 
non-relativistic limit because (i) it provides a very good numerical estimate, see Fig. 12, (ii) it allows 
us to eliminate the parametric dependence in (5. 8), (5. 10) and write f(q m ) explicitly in terms of p 
and eo, and (iii) spin effects can be accounted for [55, 56]. 

A. Arbitrary binding 

Motion of a particle can be treated non-relativistically if its momentum is much less than its 
mass. In such a case eo ~ m or ej, = m — eo <C m. Additionally, motion of a charged particle in 
electromagnetic field is non-relativistic if E <C B. Indeed, the average velocity of a non-relativistic 
particle is of order v ~ E / B = p. Thus, the non-relativistic limit is obtained by taking the limits 
Q> = £b/ m *C 1 an d p <C 1. In these limits the extremum conditions (5. 11a), (5. lib) reduce to 

t = , (5.16a) 

P 

^^ = ^. (5.16b) 
to p - q m 

Out of two solution to (5.16a) we pick the following one 

q m = -^r 2 p 2 -2e b . (5.17) 
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The sign of q m is fixed using (5.16b) by noticing that ta^ih to/to <C 1. Eliminating Q ni gives: 



Tq - (r cothr - 1) 



(5.18) 



where 



7 



(5.19) 



7 is analogous to the adiabaticity parameter of Keldysh [64]. Taking the non-relativistic limit of 
(5.12) and using (5.17) yields 

2m 2 (2e fe ) 3 / 2 



where 5(7) is the Keldysh function [64] 

9{l) 



3eE 



(5.20) 



3tq 
2 7 



7 V7 



1 



1/2' 



(5.21) 



In Fig. 12 we show the dimensionless ratio f m eE/m 2 as a function of the binding energy e?, (in 
units of m) for several values of p. The vacuum binding energy of J/ip corresponds to e& = 0.68. 
We observe an excellent agreement between the full relativistic calculation and the non-relativistic 
approximation. At p = 0.9 and = 0.68 the difference between the two lines is 10% and can be 
further improved by considering higher order corrections to f m [61]. 




FIG. 12: Dimensionless function f m eE/m 2 versus ef, for different values of p. The solid line is the full 
relativistic calculation, the dashed line is the non-relativistic approximation. J/%p binding energy in vacuum 
corresponds to = 0.68. 
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B. Weak binding 

Of special interest is the limit of weak binding 7 <S 1, i.e. y/2eb *C p. Expanding (5.18) at small 
7 and ro yields 

to = 7 (l + 7^7 2 ) (5-22) 
and substituting into (5.21) and subsequently into (5.20) yields 

fm = ~(2e b )^. (5.23) 
Hence, the quarkonium dissociation probability reads [62] 

w = exp \-3^Erj- (5 - 24) 

Since the quasi-classical approximation employed in this paper is valid if f(q m ) S> 1, it follows that 
the binding energy must satisfy 

y ! <^e b <^p 2 m. (5.25) 

Note also that we work in the approximation of the short-range binding potential meaning that 
\f2e~\j <C l/(mR), see Sec. 5 §1. 

C. Strong binding 
In the limit 7 > 1, (5.18) and (5.21) imply that 

ro = ^, 5(7) = y- (5.26) 



Substituting (5.26) into (5.20) we derive 



/ m = |§. (5.27) 



Thus, quarkonium dissociation probability in the case of strong binding is 



This formula is valid when 



p 2 m , y/eEp <C e b < 1/R 



(5.29) 



43 



D. Contribution of quark spin 

So far I have neglected the contribution of quark spin. In order to take into account the effect 
of spin interaction with the external field, we can use squared Dirac equation for a bi-spinor ip: 

[(e - eip) 2 - (p - eA) 2 - m 2 + eS • B - iea ■ E] if> = , (5.30) 

where 





a= (5.31) 



Operators T, y and a x do not commute. Therefore, in order to apply the WKB method for calcula- 
tion of the ionization probability one actually needs to square (5.30), which leads to a differential 
equation of the fourth order in derivatives. The problem becomes more tractable in the non- 
relativistic case and for crossed fields. Spin effects in crossed fields were discussed in [65] . 
With quark spin taken into account, the non-relativistic version of (5.4) becomes: 

^ [(p z + eBx) 2 + p 2 x ] - e Ex-^s-B = -e h , (5.32) 

and hence 

p 2 x = 2m {-e b + ^s-B + eEx^j - (p z + eBx) 2 , (5.33) 

where p is the quark magnetic moment and s is the projection of spin in the direction of the 
magnetic field. For a point quark, p = Hb = gfpb • The effect of quark spin on quarkonium 
dissociation probability can be taken into account by replacing £5 — > e' b = £& — ■ B in formulas 
for f m . With this replacement, all results of this section apply to a particle with spin. Note that 
effective binding energy e' h decreases if spin is parallel to the magnetic field and increases if it is 
antiparallel. In particular, in the case of weak binding 



\ - / 2 (2e b m + 2seBf' 2 
w = E ex P "3 ^neE >■ (5 - 34) 

s=±l/2 I 

Since the non-relativistic limit provide a good approximation of the full relativistic formulas, we 
will implement the quark spin dependence using the non-relativistic prescription [55, 56]. 
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§5. Effect of electric field produced in the lab frame 



A. Origin of electric field in the lab frame 



So far I have entirely neglected possible existence of electric field in the lab frame. This field, 
which we will denote by Eq, can have two origins: (i) Asymmetry of nucleon distributions in the 
colliding heavy ions, see Fig. 3(b) and (ii) Chiral Magnetic Effect (CME) [4, 66-71], which has 
recently attracted a lot of attention. In a nutshell, if a metastable P and CP-odd bubble is induced 
by axial anomaly in the hot nuclear matter, then in the presence of external magnetic field Bq the 
bubble generates an electric field which is parallel to the magnetic one. According to [67] the value 
of the electric field Eq in the bubble is 



where the sum runs over quark flavors / and it is assumed that only three lightest flavors con- 
tribute. The value of the 0-angle fluctuates from event to event. CME refers to the macroscopic 
manifestation of this effect - separation of electric charges with respect to the reaction plane. This 
effect is a possible explanation of experimentally observed charge asymmetry fluctuations [72-74]. 

No matter what is the origin of electric field in the lab frame, it averages to zero over an ensemble 
of events. We are interested to know the effect of this field on quarkonium dissociation - this is 
the problem we are turning to now [56] . 



Ionization probability of quarkonium equals its tunneling probability through the potential 
barrier. In the WKB approximation the later is given by the transmission coefficient and was 
calculated in Sec. 5 §3. In this method contribution of the quark spin can be easily taken into 
account. Another method of calculating the ionization probability, the imaginary time method 
[75, 76], was employed in [59-61]. It also yielded in the non-relativistic approximation the pre- 
exponential factor that appears due to the deviation of the quark wave function from the quasi- 
classical approximation. Such a calculation requires matching quark wave function inside and 
outside the potential barrier [62] . Extension of this approach to the relativistic case is challenging 
due to analytical difficulties of the relativistic two-body problem. Fortunately, it was argued in 
Sec. 5 §3, that the non-relativistic approximation provides a very good accuracy in the £& <C m 
region, which is relevant in the quarkiononium dissociation problem [55, 60]. 




(5.35) 



B. Quarkonium dissociation rate 
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Given the electromagnetic field in the laboratory frame -Bo, Eq, the electromagnetic field B, E 
in the comoving frame moving with velocity V is given by 



E =E | 7L (5 + Po'V x 6q) - ( 7 L - 1)^^^} 



B =B i 7 l(&o - PoV x 60) - {il ~ 1)V 



Vb c 

V 2 



(5.36a) 
(5.36b) 



where bo = Bq/Bq is a unit vector in the magnetic field direction, po = Eq/Bq = 2a|6|/37r (see 
(5.35)) and 7L = 1/Vl - V 2 . It follows from (5.36) that 



E = E ^l + 7 2 (b xV) 2 (l + Po 2 ) 
B = B ^l + 7 2 (b xV) 2 (l + p 2 ) 



(5.37a) 
(5.37b) 



Using (5.37) we find that the angle 9 between the electric and magnetic field in the comoving frame 
is 

E B 1 



cos0 = 



EB 



y/[l + 7i(bo x V) 2 (l + Po 2 )] [1 + 7 £(&o x V) 2 (l + p 2 )} 
where we used the relativistic invariance of E ■ B. 

It is useful to introduce dimensionless parameters 7, e and p as [60] 



(5.38) 



7 




1 2e b 



E 



eE / m\ 



3/2 



, , o f ( 5 - 39 ) 

m £> 

where m is quark mass and £5 is quarkonium binding energy. I will treat the quarkonium binding 
potential in the non-relativistic approximation, which provides a very good accuracy to the dissoci- 
ation rate [55, 60]. The quarkonium dissociation rate in the comoving frame in the non-relativistic 
approximation is given by [59] 



where function g reads 



3t 
2 7 



1 T 



7 V7 



a- 1 



1/2 



1 7"n 
sin ^ cos t 

37^ 



and functions P and C are given be the following formulas: 



P= 1 - 

TO 

C = exp 



(to coth t + 

70 



sinh tq cosh tq 



In + / dT 

27 7o 



TO 

7 



2) sin 2 + sinh 2 r cos 2 



I-1/2 



£(t) t - t 



1 



(Tq T 2 ) 2 COS^ + Tq 



/ cosh To — COS T 
\ sinh tq 



/ sinh t t V 
V sinh t T / 



1/2 



sin 



(5.40) 
(5.41) 

(5.42) 
(5.43) 

(5.44) 
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The contribution of quark spin is taken into account by replacing e& — >■ e' b = e b — ■ B [55]. 
Function g represents the leading quasi-classical exponent, P is the pre- factor for the S- wave state 
of quarkonium and C accounts for the Coulomb interaction between the valence quarks. Parameter 
ro satisfies the following equation 

r 2 - sin 2 6{t cothr - l) 2 = j 2 (5.45) 

which establishes its dependence on 9 and 7. Note, that in the limit E — > the dissociation rate 
(5.40) exponentially vanishes. This is because pure magnetic field cannot force a charge to tunnel 
through a potential barrier. 

In the case that mechanism (i) is responsible for generation of electric field, Eq is the field 
permitting the entire plasma in a single event. Event average is then obtained by averaging 
(5.40) over an ensemble of events. In the case that mechanism (ii) is operative, averaging is more 
complicated. Eq. (5.40) gives the quarkonium dissociation rate in a bubble with a given value of 
0. Its derivation assumes that the dissociation process happens entirely inside a bubble and that 
is constant inside the bubble. Since in a relativistic heavy ion collision many bubbles can be 
produced with a certain distribution of 0's (with average (G) = 0) more than one bubble can affect 
the dissociation process. This will result in a distractive interference leading to reduction of the 
CP-odd effect on quarkonium dissociation. However, if a typical bubble size Rq is much larger than 
the size of quarkonium Rj, then the dissociation is affected by one bubble at a time independently 
of others, and hence the interference effect can be neglected. In this case (2.27) provides, upon a 
proper average, a reasonable estimate of quarkonium dissociation in a heavy ion collision. We can 
estimate the bubble size as the size of the sphaleron, which is of the order of the chromo-magnetic 
screening length ~ l/g 2 T, whereas the quarkonium size is of the order a s /e&. Consequently, at 
small coupling and below the zero-field dissociation temperature (i.e. when ej, is not too small) 
Ro is parametrically much larger than Rj. A more quantitative estimate of the sphaleron size is 
R ~ 1.2/a s N c T ~ 0.4 fm [77]; whereas for J/i/> Rj ~ a s /e b ~ 0.1 - 0.2 fm. Thus, based on this 
estimate bubble interference can be neglected in the first approximation. However, since the ratio 
Rj/Ro is actually not so small this effect nevertheless warrants further investigation. 

To obtain the experimentally observed J/ip dissociation rate we need to average (2.27) over the 
bubbles produced in a given event and then over all events. To this end it is important to note 
that because the dissociation rate depends only on p\ it is insensitive to the sign of the Eq field 
or, in other words, it depends only on absolute value of © but not on its sign. Therefore, it stands 
to reason that although the precise distribution of 0's is not known, (2.27) gives an approximate 
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event average with parameter representing a characteristic absolute value of the theta- angle. 



Before I proceed with the numerical calculations, let us consider for illustration several limiting 
cases. If quarkonium moves with non-relativistic velocity, then in the comoving frame electric and 
magnetic fields are approximately parallel 9 « 0, whereas in the ultra-relativistic case they are 
orthogonal 6 ~ tt/2, see (5.38). In the later case the electromagnetic field in the comoving frame 
does not depend on E$ as seen in (5.37) and therefore the dissociation rate becomes insensitive 
to the CME. In our estimates I will assume that po < 1 which is the relevant phenomenological 
situation. Indeed, it was proposed in [67] that po ~ a <C 1 produces charge fluctuations with 
respect to the reaction plane of the magnitude consistent with experimental data. 

1) ^ 0, i.e. electric and magnetic fields are approximately parallel. This situation is realized 
in the following two cases, (i) Non-relativistic quarkonium velocities: V <C po or (ii) motion of 
quarkonium at small angle <j> to the direction of the magnetic field &o : <fi ^ Po/jlV- In both cases 
E ~ Eq and B ~ Bq. This is precisely the case where the dissociation rate exhibits its strongest 
sensitivity to the strength of the electric field Eq generated by the local parity violating QCD 
effects. Depending on the value of the 7 parameter defined in (5.39) we can distinguish the case 
of strong electric field 7> 1 and weak electric field 7 <ti 1 [61]. In the former case, g = (3/8)7, 
P = (8/e) 1 / 2 7e-T 2 / 2 , C = e^/ 2 /^. Substituting into (2.27) the dissociation rate reads 



where the electromagnetic field in the comoving frame equals one in the laboratory frame as was 
mentioned before. 



C. Limiting cases 




(5.46) 



In the later case, g = P = C = 1 and 




(5.47) 



2) ~ vr/2, i.e. electric and magnetic fields are approximately orthogonal."'' This occurs for an 



ultra-relativistic motion of quarkonium V — > 1. In this case 




(5.48) 



* Note, that the limit 7 3> 1 is different in = | and 9 < | cases [59]. 
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This case was discussed in detail in our previous paper [55]. In particular for 7 <C 1 we get 

(5.49) 

eE \m J 

Due to (5.46) and (5.49) dependence of w on Eq is weak unless po S> 1. 



§6. Dissociation rate of J/ip 

One of the most interesting applications of this formalism is calculation of the dissociation rate 
of J/ if) which is considered a litmus test of the quark-gluon plasma [54]. Let z be the heavy ions 
collision axis; heavy-ion collision geometry implies that b^ - z = 0. The plane containing z-axis and 
perpendicular to the magnetic field direction is the reaction plane. We have 



(bo x V) 2 = V 2 + Vl sin 2 



(5.50) 



where 4> is the angle between the directions of Bq and V± and I denoted vector components in the 
xy-plane by the subscript _L. We can express the components of the quarkonium velocity V in 
terms of the rapidity rj as V z = tanhry, V± = P±/(M± cosh 77), where P and M are the quarkonium 
momentum and mass and M 2 = M 2 + P 2 . 
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FIG. 13: Dissociation rate of J/ip at eBo = 15mJ, (f> = 7r /2 (in the reaction plane), 77 = (midrapidity) as 
a function of (a) P± at St = 0.16 GeV and (b) at P± = 1 GeV. 

Results of numerical calculations are exhibited in Figs. 13-15 [56]. In Fig. 13 I show the 
dissociation rate of J/tp for several values of the electric field Eq induced by the Chiral Magnetic 
Effect. Note, that the typical size of the medium traversed by a quarkonium in magnetic field 
can be estimated very conservatively as a few fm. Therefore, w ~ 0.3 — 0.5 fm _1 corresponds to 



complete destruction of J/ip's. This means that in the magnetic field of strength eB$ ~ 15m 2 
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all J/ip's with P± > 0.5 GeV are destroyed independently of the strength of Eq. Since magnetic 
field strength decreases towards the QGP periphery, most of J/ip surviving at later times originate 
from that region. Effect of electric field Eq is strongest at low P±, which is consistent with our 
discussion in the previous section. The dissociation rate at low P± exponentially decreases with 
decrease of Eq. Probability of quarkonium ionization by the fields below Eq < 0.1-Bo (i- e - Po ^5 0.1) 
is exponentially small. This is an order of magnitude higher than the estimate po ~ a of electric 
field due to CME effect as proposed in [67]. 




FIG. 14: Contour plot of the dissociation rate of J/ip as a function of Eb and eEg at eBo = 15m^, <p = ir/2 
(in the reaction plane), 77 = (midrapidity) and P± =0.1 GeV. Numbers inside boxes indicate the values 
of w in 1/fm. 



As the plasma temperature varies, so is the binding energy of quarkonium although the precise 
form of the function £&(T) is model-dependent. The dissociation rate picks at some e® < e^ ac (see 
Fig. 13(b)), where e^ ac is the binding energy in vacuum, indicating that J/ip breaks down even 
before Eb drops to zero, which is the case at Bq = 0. This is a strong function of Eq as can be 
seen in Fig. 14. It satisfies the equation dw/deb = 0. In the case 7^1 (5.47) and (5.49) imply 
that 

e ° = tM ■ ^ (551) 

At 7 » 1 and 9 = tt/2 we employ (5.46) to derive the condition (e°) 2 + eBe\j2m — eE 2 /B = 0. 
In view of (5.48) E ~ B and we obtain 



£ b 




'eB. 



7 > 1 



(5.52) 
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where in the last step I used that eB <C m 2 . For a given function £&(T) one can convert into 
the dissociation temperature, which is an important phenomenological parameter. 
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FIG. 15: (a) Angular distribution of J/tp dissociation rate at eBg — 15mJ, rj = at different Eq and Pj_ 
(in GeV's). Magnetic field Bo points in the positive vertical direction. Reaction plane coincides with the 
horizontal plane, (b) Rescaled second Fourier-harmonic vi of the azimuthal distribution as a function of 
P±_ . (P) is the azimuthal average of the survival probability and t is the time spent by J ftp in the P-odd 
bubble. 



In the absence of electric fied Eq, the dissociation probability peaks in the direction perpendic- 
ular to the direction of magnetic field bo, i.e. in the reaction plane. Dissociation rate vanishes in 
the &o direction. Indeed, for V ■ &o = (5.37) implies that E = 0. This feature is seen in the left 
panel of Fig. 15. At finite Eq the dissociation probability is finite in the bo direction making the 
azimuthal distribution more symmetric. The shape of the azimuthal distribution strongly depends 
on quarkonium velocity: while at low V the strongest dissociation is in the direction of the reaction 
plane, at higher V the maximum shifts towards small angles around the &o direction. Extrema of 
the azimuthal distribution are roots of the equation dw/dcj) = 0. At 7 <C 1 it yields minimum at 
00 = 0, maximum at 4>q = ir/2 and another maximum that satisfies the condition (neglecting the 
spin-dependence of 

eEo^l + 7 £(K 2 + VI sin 2 + p^) = (5.53) 

In order to satisfy (5.53) (fio must decrease when V increases and visa versa. This features are seen 
in the left panel of Fig. 15. 

Spectrum of quarkonia surviving in the electromagnetic field is proportional to the survival 
probability P = 1 — wt, where t is the time spent by the quarkonium in the field. Consider P as 
a function of the angle \ between the quarkonium velocity and the reaction plane x = 71 / 2 — (j). 
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Fourier expansion of P in \ reads 



P(X) = -P + ^PnCOs(n X ), P, 



n=l 



n 



7T 



1 



/_ 



— 7T 



■7T 



P(x) cos(nx) c?x • 



(5.54) 



Ellipticity of the distribution is characterized by the "elliptic flow" coefficient v 2 defined as 



where (P) denotes average of P over the azimuthal angle. These formulas are applicable only as 
long as wt < 1 because otherwise there are no surviving quarkonia. In the right panel of Fig. 15 
[56] I show V2 (P) /t, which is independent of t, as a function of P±. As expected, in the absence of 
the CME, t>2 is negative at low P± and positive at high P±. V2 changes sign at P± that depends on 
the strength of the electric field. It decreases as Eq increases until at Eq — -Bo it becomes positive 
at all P±. Fig. 15(b) provides the low bound for V2 because (P) < 1 and t > 1 fm . We thus expect 
that magnetic field strongly modifies the azimuthal distribution of the produced J/ip's. Role of the 
magnetic field in generation of azimuthal anisotropies in heavy-ion collisions has been pointed out 
before in [27, 52]. 

In summary, we observed that J/ip dissociation energy increases with magnetic field strength 
and quarkonium momentum. As a consequence, quarkonia dissociate at lower temperature than 
one would have expected based on calculations neglecting magnetic field [55, 56]. Fig. 13 indicates 
that in heavy-ion collisions at the LHC, all J/ip's moving with P± > 0.5 GeV in the reaction plane 
would dissociate with probability of order unity even if the QGP effect were completely negligible. 
If electric field fluctuations shown in Fig. 3 are taken into account, then even low P± J/i/j's are 
destroyed. However, Chiral Magnetic Effect has negligible effect on J /if} dissociation. 

Although magnetic fields in pp and pA collisions are much weaker than in AA collisions, they are 
still strong enough to cause J/tp dissociation at sufficiently high momenta P±. A truly spectacular 
feature of such process would be J/ip decay into two heavier D-mesons. 

The effect of J/tp dissociation in a magnetic field vanishes in the direction parallel to the 
magnetic field, i.e. perpendicular to the reaction plane. Therefore, J/ip dissociation gives negative 
contribution to the total azimuthal asymmetry coefficient V2- It remarkable that presence of electric 
field reverses this effect making V2 positive. 




(5.55) 
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6. ELECTROMAGNETIC RADIATION BY QUARK-GLUON PLASMA IN MAGNETIC 

FIELD 

§1. Necessity to quantize fermion motion 

In Sec. 3 we discussed synchrotron radiation of gluons by fast quarks. Our main interest was 
the energy loss problem. In this section we turn to the problem of electromagnetic radiation by 
QGP, viz. radiation of photons by thermal fermions [78]. In this case quasi- classical approximation 
that we employed in Sec. 3 and Sec. 4 is no longer applicable and one has to take into account 
quantization of fermion motion in magnetic field. 

Electromagnetic radiation by quarks and antiquarks of QGP moving in external magnetic field 
originates from two sources: (i) synchrotron radiation and (ii) quark and antiquark annihilation. 
QGP is transparent to the emitted electromagnetic radiation because its absorption coefficient is 
suppressed by a 2 . Thus, QGP is shinning in magnetic field. The main goal of this paper is to 
calculate the spectrum and angular distribution of this radiation. In strong magnetic field it is 
essential to account for quantization of fermion spectra. Indeed, spacing between the Landau levels 
is of the order eB/e (e being quark energy), while their thermal width is of the order T. Spectrum 
quantization is negligible only if eB/e <C T which is barely the case at RHIC and certainly not the 
case at LHC (at least during the first few fm's of the evolution). Fermion spectrum quantization is 
important not only for hard and electromagnetic probes but also for the bulk properties of QGP. 

§2. Synchrotron radiation 

Motion of charged fermions in external magnetic field, which I will approximately treat as 
spatially homogeneous, is quasi-classical in the field direction and quantized in the reaction plane, 
which is perpendicular to the magnetic field and span by the impact parameter and the heavy 
ion collision axis. In high energy physics one usually distinguishes the transverse plane, which 
is perpendicular to the collision axis and span by the magnetic field and the impact parameter. 
In this section I use notation in which three-vectors are discriminated by the bold face and their 
component along the field direction by the plain face. Momentum projections onto the transverse 
plane are denoted by subscript _L. 

In the configuration space, charged fermions move along spiral trajectories with the symmetry 
axis aligned with the field direction. Synchrotron radiation is a process of photon 7 radiation by a 
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fermion / with electric charge e/ = Zfe in external magnetic field B: 



f(e f ,j,p) -»• f(e f ,k,q) +j(k) 



(6.1) 



where k is the photon momentum, p, q are the momentum components along the magnetic field 
direction and indicies j, k = 0, 1, 2, . . . label the discrete Landau levels in the reaction plane. The 
Landau levels are given by 



In the constant magnetic field only momentum component along the field direction is conserved. 
Thus, the conservation laws for synchrotron radiation read 



where uj is the photon energy and 9 is the photon emission angle with respect to the magnetic field. 
Intensity of the synchrotron radiation was derived in [79] . In [80-83] it was thoroughly investigated 
as a possible mechanism for 7-ray bursts. In particular, synchrotron radiation in electromagnetic 
plasmas was calculated. Spectral intensity of angular distribution of synchrotron radiation by a 
fermion in the j'th Landau state is given by 



where Tj k = (1 + Sjo)(l + 5 k o) accounts for the double degeneration of all Landau levels except the 
ground one. The squares of matrix elements M., which appear in (6.4), corresponding to photon 
polarization perpendicular and parallel to the magnetic field are given by, respectively, 




(6.2) 



£j = uj + e k , p = q + uj cos 9 , 



(6.3) 




(6.4) 



4eje k \M± 



(ejEk-pq - m 2 )[ij fc _! + /|_ ljfe ] + 2 y /2je f By/2ke f B[I jtk _ 1 I j _ 1>k ] . 



(6.5) 



cos 2 9{(s j e k -pq-m 2 )[ll k _ 1 + l]_ ljk ] - 2^2j^B^2ke f B[I j ^ 1 I j _ 1 , k }} 



- 2cos9sm9{py / 2ke f B[I j - ljk I j ^ ltk - 1 + I jtk - 1 I j)k ] 



+ q^2je f B[I j 

,klj-l,k + Ij-l,k-llj,k-l]} 



+ sin 2 9{(£ j£k +pq- m 2 )^^! + / 2 fc ] + 2^2je f B^/2ke f B{I j _ l ^ l I h k )} 



(6.6) 



where for j > k, 




(6.7) 
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and Ij,k{x) = when k > j. (Ij-i are identically zero). The functions L J k ~ h (x) are the 

generalized Laguerre polynomials. Their argument is 

x = sin 2 9 . (6.8) 

2e f B { ' 

Angular distribution of radiation is obtained by integrating over the photon energies and re- 
membering that Ek also depends on to by virtue of (6.2) and (6.3): 

/ fe=0 J 

where photon energy uj is fixed to be 

u>* = {( £j -pcos9) - [( £j -pcos9) 2 - 2e f B(j - A;) sin 2 0] 1/2 } . (6.10) 

In the context of heavy-ion collisions the relevant observable is the differential photon spectrum. 
For ideal plasma in equilibrium each quark flavor gives the following contribution to the photon 
spectrum: 

diV synch roc B(2N C )V ^ ^ dP r w, m , , 

/ J 00 j=0fc=0 

where 2iV c accounts for quarks and antiquarks each of N c possible colors, and (2 — 5jfl) sums over 
the initial quark spin. Index / indicates different quark flavors. V stands for the plasma volume. 
The statistical factor /(e) is 

m-^i- (•■«) 

The 5-function appearing in (6.4) can be re-written using (6.2) and (6.3) as 



where 



p*± =| cos9(rrij — m k + oj sin 9) 

± \J [{rrij + m,&) 2 — oj 2 sin 2 — m^) 2 — oj 2 sin 2 0] j/(2wsin 2 #). (6.14) 
The following convenient notation was introduced: 

m 2 = m 2 + 2je f B, m\ = m 2 + 2ke f B . (6.15) 
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The physical meaning of (6.14) is that synchrotron radiation of a photon with energy u> at angle 
9 by a fermion undergoing transition from j'th to fc'th Landau level is possible only if the initial 
quark momentum along the field direction equals p± . 

Another consequence of the conservation laws (6.3) is that for a given j and k the photon energy 
cannot exceed a certain maximal value that will be denoted by uj s ,jk- Indeed, inspection of (6.14) 
reveals that this equation has a real solution only in two cases 

(i) rrij — rrik > uj sin 9 , or (ii) rrij + < uj sin 9 . (6.16) 

The first case is relevant for the synchrotron radiation while the second one for the one-photon pair 
annihilation as discussed in the next section. Accordingly, allowed photon energies in the j — > k 
transition satisfy 

rrij - m k ^Jm 2 + 2je f B - ^m 2 + 2ke f B 
u < ^s,jk = . Q = ~ r-^ • (6.17) 

No synchrotron radiation is possible for co > uj s jh- I n particular, when j = k, uj s jk = 0, i.e. no 
photon is emitted, which is also evident in (6.10). The reason is clearly seen in the frame where 
p = 0: since £j > constraints (6.2) and (6.3) hold only if oj = 0. 

Substitution of (6.4) into (6.11) yields the spectral distribution of the synchrotron radiation 
rate per unit volume 

dArsynch 2N c zja r S(p-p* ± ) 

/ j=0 fc=0 J ± \ej e k \ 

x{\M ± \ 2 + \M {l \ 2 } f(ej)[l -f(e k )}, (6.18) 

where # is the step-function. 

The natural variables to study the synchrotron radiation are the photon energy ui and its emis- 
sion angle 9 with respect to the magnetic field. However, in high energy physics particle spectra are 
traditionally presented in terms of rapidity y (which for photons is equivalent to pseudo-rapidity) 
and transverse momentum k±. k± is a projection of three- momentum k onto the transverse plane. 
These variables are not convenient to study electromagnetic processes in external magnetic field. 
In particular, they conceal the azimuthal symmetry with respect to the magnetic field direction. 
To change variables, let z be the collision axis and y be the direction of the magnetic field. In 
spherical coordinates photon momentum is given by k = u(sm a cos <px + sin a sin (fiy + cosai), 
where a and 4> are the polar and azimuthal angles with respect to z-axis. The plane xz is the 
reaction plane. By definition, k ■ y = cos 6 implying that cos 9 = sin a sin (p. Thus, 

fTn cos 9 , a , 

k± = Jkl + k 2 = —— , y = - In tan - . 6.19 
V y sm0 2 
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The second of these equations is the definition of (pseudo)-rapidity. Inverting (6.19) yields 

sin (b 

uj = k± cosh y , cos 9 = — . (6.20) 

coshy 

Because dy = dk z /uj the photon multiplicity in a unit volume per unit time reads 



dAr synch 



diY synch 



dAr synch 



UJ 



dVdt d 2 k± dy dVdt d 3 k dVdt uduodQ. 



(6.21) 



0.001 



0.001 s 




k T . GeV 



(a) 



(b) 



FIG. 16: Spectrum of synchrotron radiation by u quarks at eB = m^, y = 0, <j> = 7r/3: (a) contribution of 
10 lowest Landau levels j < 10; several cutoff frequencies are indicated; (b) summed over all Landau levels. 
ni u = 3 MeV, T = 200 MeV. Adopted from [78]. 

Fig. 16 displays the spectrum of synchrotron radiation by u quarks as a function of k± at fixed 
(p [78]. At midrapidity y = (6.20) implies that k± = uj. Contribution of d and s quarks is 
qualitatively similar. At eB 3> m 2 , quark masses do not affect the spectrum much. The main 
difference stems from the difference in electric charge. In panel (a) only the contributions of the 
first ten Landau levels are displayed. The cutoff frequencies u} s ,jk can be clearly seen and some of 
them are indicated on the plot for convenience. The azimuthal distribution is shown in Fig. 17. 
Note, that at midrapidity <j) = tt/2 — 9. Therefore, the figure indicates that photon production 
in the direction of magnetic field (at = 7r/2) is suppressed. More photons are produced in the 
direction of the reaction plane = 0. This results in the ellipticity of the photon spectrum that 
translates into the positive "elliptic flow" coefficient V2- It should be noted, that the classical 
synchrotron radiation has a similar angular distribution. 

In order to compare the photon spectrum produced by synchrotron radiation to the photon 
spectrum measured in heavy-ion collisions, the u, d and s quarks contributions were summed up. 
Furthermore, the experimental data from [84] was divided by Vt, where t is the magnetic field 
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FIG. 17: Azimuthal distribution of synchrotron radiation by it-quarks at k± = 0.2 GeV, eB = m^, y = 0. 
m u = 3 MeV. Adopted from [78] 

relaxation time. The volume of the plasma can be estimated as V = irR 2 t with R ~ 5 fm being 
the nuclear radius. Therefore, 

gVgcp _ rfiVclp 1 _ diVclp /GeV\ 4 / lfm 
dVdtd 2 k ± dy d 2 k ± dy ttRH 2 d 2 k ± dy \ 14.9 / V t 

The results are plotted in Fig. 18. In panel (a) it is seen that synchrotron radiation gives a 
significant contribution to the photon production in heavy-ion collisions at RHIC energy. This 
contribution is larger at small transverse momenta. This may explain enhancement of photon 
production observed in [84]. Panel (b) indicates the increase of the photon spectrum produced by 
the synchrotron radiation mechanism at the LHC energy. This increase is due to enhancement of 
the magnetic field strength, but mostly because of increase of plasma temperature. This qualitative 
features can be better understood by considering the limiting cases of low and high photon energies. 

One possible way to ascertain the contribution of electromagnetic radiation in external magnetic 
field is to isolate the azimuthally symmetric component with respect to the direction of the magnetic 
field. It seems that synchrotron radiation dominates the photon spectrum at low k±. Thus, 
azimuthal symmetry can be easily checked by simply plotting the multiplicity vs uj, 9 and tp, where 
uj is photon energy, 6 is emission angle with respect to the magnetic field and tp is azimuthal angle 
around the magnetic field direction (which is perpendicular both to the collision axis and to the 
impact parameter). In Fig. 16(a) it is also seen that in these variables it may be possible to discern 
the cutoff frequencies uj s jk that appear as resonances (in Fig. 16 y = so k± = uj). Note that 
averaging over the azimuthal angle a around the collision axis direction destroys these features, 
see Fig. 18. 




(6.22) 
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FIG. 18: Azimuthal average of the synchrotron radiation spectrum of u,d,s quarks and their corresponding 
antiquarks. (a) eB = m^, y — compared to the experimental data from [84] divided by Vt = 25-7T fm 4 
(dots) and Vt = 9x 25tt fm 4 (stars), (b) eB = m£, T = 200 MeV, y = (solid line) compared to eB = \hm\ , 
T = 400 MeV, y = (dashed line). m u = 3 MeV, m d = 5 MeV, m s = 92 MeV. Adopted from [78]. 

A. Low photon energy 



The low energy part of the photon spectrum satisfies the condition ui -C ^/ejB. The corre- 
sponding initial quark momentum component along the field p and energy 8j follow from (6.14) 
and (6.2) and are given by 

(j - k)e f B(cos6±l) 



P± 



oj sin 



+ 0(u), Ej ^\p* ± \ + 0(<j). 



(6.23) 



Evidently, e,- 3> e5. In practice, magnetic field strength satisfies yeB > T, so that ej 3> T. 
Therefore, synchrotron radiation is dominated by fermion transitions from low Landau levels due 
to the statistical factors appearing in (6.11). 

For a qualitative discussion it is sufficient to consider the 1 — > transition. In this case the 
matrix elements (6.5) and (6.6) read 



\M 



1,0|2 



— \ if (eie - pq cos 2 9 - m 2 ) + cos 8 sin 9q^/2ejBI lfi Io,o\ ■ (6.24) 
i£n l ' > 



2eie 



Assuming that the field strength is supercritical, i.e. efB 3> m?, but keeping all powers of u (for 
future reference) (6.14) reduces to 



p* ± « ^ {2e f £(cos 6> ±l)+w 2 sin 2 (9(cos 19 =Fl)} . 

2w sin 9 



(6.25) 
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Furthermore, using the conservation laws (6.3) we obtain in this approximation 

1 



|2e/.B(cos6»± 1) - w 2 sin 2 <9(cos 6> =F 1)| , 



2oj sin 2 9 

q± = ^-^(2e f B - uj 2 sin 2 6) (cos 6 ± 1) , 

2u> sin 9 

£o± =\q\ ■ 

The values of the non- vanishing matrix elements I~ \. defined by (6.7) are 



(6.26) 
(6.27) 
(6.28) 



h,o(x) 



-x l ' 2 e- x ' 2 



Io,o(x) 



-x/2 



For j = 1, k = we write using (6.17) w S) io = y/^ejB/ sin 9. Then (6.8) implies x 
Substituting (6.25)-(6.29) into (6.24) gives 
1 



u 2 /u> 2 



(6.29) 

's,10- 



1,0,2 
± I 



-xe 



cos6>(l + x)±(l-x) ,„ 2(1 -x) cos 6> sin 2 (9 
1 tttz 4—4: r cos^ 9 - 



(6.30) 



cos0(l-x)±(l + x) cos0(l - x) ± (1 + x)_ 

According to (6.18) the contribution of the 1—^0 transition to the synchrotron radiation reads 
[78] 

^ysynch,10 „ - 2 \ r z'i<\ 

M _ ms ft 4- (\ 4- r\ 

(6.31) 



VdtdVLduj 



,io 2iV c z?a e f B ^ , n 

E7 -£^ F ^£iE/ta)[i-/fo.)llA# , l 

f ± 



(l-x)cos6>±(l + x) . , 



-2x(cos6>T 1) 

Consider radiation spectrum at = 7r/2, i.e. perpendicular to the magnetic field. The spectrum 
increases with x and reaches maximum at x = 1. Since x = uj 2 /(2efB), spectrum decreases with 
increase of B at fixed uj. This feature holds at low x part of the spectrum for other emission angles 
and even for transitions form higher excited states. However, at high energies, it is no longer 
possible to approximate the spectrum by the contribution of a few low Landau levels. In that 
case the typical values of quantum numbers are j,k 3> 1. For example, to achieve the numerical 
accuracy of 5%, sum over j must run up to a certain j max . Some values of j max are listed in Table I 
[78]. 



B. High photon energy 

The high energy tail of the photon spectrum is quasi-classical and approximately continuous. In 
this case the Laguerre polynomials can be approximated by the Airy functions or the corresponding 
modified Bessel functions. The angular distribution of the spectrum can be found in [80]: 



^ synch _ z _h n f um 2 j e f BT sm9 /T 

Vdtdndio 4- vr 4T3 V m 3 ' [ ' 
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30 
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8 12 16 



TABLE I: The upper summation limit in (6.18) that yields the 5% accuracy. j max is the highest Landau 
level of the initial quark that is taken into account at this accuracy. Throughout the table y = 0. 

provided that uj 3> mJ mT/efB sin 0. Here n/ is number density of flavor /, which is independent 
of B: 

n^'-^tf^^TK (,33) 

Here summation over j was replaced by integration. It follows that this part of the spectrum 
increases with magnetic field strength as \[B and and with temperature as \/Te~ UJ / T . Therefore, 
variation of the spectrum with T is much stronger than with B. The T dependence is shown in 
Fig. 19. 



0.001 - 




k T , GeV 

FIG. 19: Variation of the synchrotron spectrum with plasma temperature. Lower line: T — 200 MeV, upper 
line: T = 250 MeV. Other parameters are the same as in Fig. 18(a). Adopted from [78]. 

Unlike time-dependence of magnetic field, time-dependence of temperature is non-negligible even 
during the first few fm/c. Final synchrotron spectrum, which is an average over all temperatures, 
is dominated by high temperatures/early times. However, the precise form of time-dependence of 
temperature is model-dependent. Therefore, the spectrum is presented at fixed temperatures, so 
that a reader can appreciate its qualitative features in a model-independent way. 
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§3. Pair annihilation 

The theory of one-photon pair annihilation was developed in [85, 86]. It was shown in [87] that 
in the super-critical regime eB 3> m 2 one-photon annihilations is much larger than the two-photon 
annihilation. In this section the one-photon annihilation of q and q pairs in the QGP is calculated. 

For qq pair annihilation the conservation of energy and momentum is given by 

£j+£fc=w, p + q = ucosd. (6.34) 

The spectral density of the annihilation rate per unit volume reads 

^ annih ^ azfcNc ™™ f 2e f B f 2e f B 

/ J j=0k=0 J J 

xS{p + q - oj cos 9)5{e j + e k - uj){\T±\ 2 + |7jj | 2 } , (6.35) 

where the matrix elements T can be obtained from (6. 5), (6. 6) by making substitutions e k — > —e k , 
q — > —q and are given by 

4^171 1 2 =(e 3 e k -pq + m 2 )[ij fc _! + / 2 _ 1)fe ] - 2 y /2j^B y /2ke f B[I jtk _ 1 I j _ 1>k ] . (6.36) 
4e i£fc |^| 2 =cos 2 ^{(e jefe - M + m 2 )[/ 2 fc _! +/ 2 _ ljfc ] + 2y/2je f B y /2ke f B[I jtk - 1 I j - 1>k ]} 
- 2cos6>sin0{ - py/2kefB[Ij- lik Ij-i,k-i + Ij,k-ilj,k] 



+ qy / 2je f B[I j)k I j -i )k + / ? --i,fc-ii J -,A : -i] } 

+ sin 2 0{( £j£fc +pg + m 2 )^^! + I 2 *] - 2 y /2je f B y /2ke f B(I j - hk - 1 I jtk )} , 

(6.37) 

with the same functions Jjj as in (6.7). Integration over g removes the delta function responsible 
for the conservation of momentum along the field direction. The remaining delta function is 
responsible for energy conservation and can be written in exactly the same form as in (6.13) with 
particle energies and momenta now obeying the conservation laws (6.34). It is straightforward to see 
that momentum p± is still given by (6. 14), (6. 15). The photon spectrum produced by annihilation 
of quark in state j with antiquark in state k has a threshold uj a ^j that is given by the case (ii) in 
(6.16): 

mj + m k ^m 2 + 2je f B + ^m 2 + 2ke f B 

U > W a ij — — — — — . (6.38J 

sin 6 sin 

Thus, the spectral density of the annihilation rate per unit volume is 

djV annih azjioN c ^^ Q , ,[ ^ 5(p - p* ± ) 



VdKlvm E '/"EE" U / ^E|Z_Xi 

/ j=0 k=0 J ± I £j £k 



X 



{\n\ 2 + \T\\\ 2 }f(e 3 )f(e k ). (6.39) 
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Passing to y and p± variables in place of u and is similar to (6.21). 

The results of the numerical calculations are represented in Fig. 20. Panel (a) shows the spec- 
trum of photons radiated in annihilation of u and u. We conclude that contribution of the annihi- 
lation channel is negligible as compared to the synchrotron radiation. 




Pt, CeV 

(a) (b) 

FIG. 20: Photon spectrum in one-photon annihilation of u and u quarks. eB — mj, y — 0. (a) fcj^-spectrum 
at <j> = 7r/3, (b) azimuthal angule distribution at k± = 1 GeV. Adopted from [78]. 

In summary, results of the calculations presented in this section indicate that photon produc- 
tion by QGP due to its interaction with external magnetic field give a considerable contribution 
to the total photon multiplicity in heavy-ion collisions. This is seen in Fig. 18 were the model 
calculation is compared with the experimental data [84]. The largest contribution to the photon 
multiplicity arises from photon momenta of the order of V eB. This may provide an explanation 
of the photon excess observed by the PHENIX experiment [84]. Similar mechanism may also be 
responsible for enhancement of low mass di-lepton production that proceeds via emission of virtual 
photon which subsequently decays into di-lepton pair. 



63 



7. SUMMARY 

Analytical and numerical calculations indicate existence of extremely powerful electromagnetic 
fields in relativistic heavy ion collisions. They are the strongest electromagnetic fields that exist in 
nature. They evolve slowly on characteristic QGP time scale and therefore have a profound effect 
on dynamics of QGP. In this review I described the recent progress in understanding of particle 
production in presence of these fields. Treating the fields as quasi-static and spatially homogeneous 
allowed us to use analytical results derived over the past half-century. This is however the main 
source of uncertainty that can be clarified only in comprehensive numerical approach based on 
relativistic magnetohydrodynamics. 

I discussed many spectacular effects caused by magnetic field. All of them have direct phe- 
nomenological relevance. Breaking of spherical symmetry by magnetic field in the direction per- 
pendicular to the collision axis results in azimuthal asymmetry of particle production in the reac- 
tion plane. Fast quarks moving in magnetic field radiate a significant fraction of their energy. All 
electromagnetic probes are also naturally affected by magnetic field. Therefore, all experimental 
processes that are being used to study the properties of QGP have strong magnetic field depen- 
dence. In addition, the QCD phase diagram is modified by magnetic field as has been extensively 
studied using model calculations [88-109] and lattice simulations [110-122]. Entanglement of ef- 
fects produced by magnetic field with conventional QGP ones makes it difficult to quantify the role 
of magnetic field in QGP dynamics. A unique observable is polarization of leptons escaping from 
QGP, which can be induced only by magnetic field, see Sec. 3 §3. 

Profound influence of magnetic field on properties of QGP is truly remarkable. Hopefully, 
progress in theory will soon be matched by experimental investigations that will eventually discover 
properties of QCD at high temperatures and strong electromagnetic fields. 
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